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ABSTRACT 

A  theoretical  study  and  computer  analysis  of  the  cathode 
spot  of  a  unipolar  arc  was  conducted.  The  underlying  theories 
of  Plasma  Physics,  Space  Charge  Effects,  and  Electron  Emission 
needed  for  an  understanding  of  cathode  phenomena  are  presented. 
Two  models  of  unipolar  arcing  are  reviewed  and  an  analysis  of 
the  cathode  spot  is  begun.  A  stationary  model  of  the  cathode 
spot  is  formulated  with  a  system  of  equations  that  is  not 
closed.  A  method  of  solution,  using  the  Steenbeck  Minimum 
Principle,  is  utilized  and  a  computer  program  developed  to 
determine  the  cathode  spot  parameters.  Results  of  the  analysis 
are  presented  for  arc  currents  of  100,  150,  175,  200,  300,  and 
400  Amperes.  It  is  found  that  for  copper  stable  arcs  occur 
only   for   arc   currents    greater   than  200    Amperes. 
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I.    INTRODPCTTON 

The  study  of  cathode  phenomena  in  gaseous  discharges  and 
electrical  arcs  has  progressed  over  the  last  hundred  years.  In 
the  last  twenty  years,  this  research  has  been  revived  as  the 
interest  in  high  power  devices  for  pulsed  power  applications, 
switching  devices,  and  nuclear  fusion  power  generation  has 
increased.  The  discovery  in  1958  of  a  new  form  of  electrical 
arc,  one  that  burns  between  a  metal  plate  and  a  dense  plasma, 
called  a  unipolar  arc  has  further  increased  interest  in  the 
cathode  phenomena.  For  both  the  bipolar  and  unipolar  arcs,  the 
current  transfer  occurs  in  localized  regions  between  a  cathode 
and  plasma  which  is  referred  to  as  the  cathode  spot.  It  has 
become  increasingly  apparent  that  the  phenomena  occurring  in 
the  cathode  spot  of  a  unipolar  arc  and  that  occurring  in  a 
bipolar  arc  are  similar.  Therefore,  much  of  the  theoretical  work 
done  on  the  cathode  spot  of  a  bipolar  arc  can  be  applied 
towards  the  understanding  and  formulation  of  a  consistent 
theory   of   the    unipolar   arc. 

This  study  will  begin  with  a  short  review  of  some  basic 
theories  of  plasma  physics,  space  charge  effects,  and  electron 
emission  mechanisms.  Two  useful  unipolar  arcing  models  are  then 
reviewed.  The  stationary  cathode  spot  will  then  be  studied  and 
a    set   of   equations    presented    that   can   be    used   to    determine    it's 


parameters.  It  will  be  found  that  the  system  of  equations  is 
not  closed,  however,  a  method  for  the  solution  of  the  system 
can  be  formulated  by  using  the  Steenbeck  Minimum  Principle. 
Finally,  results  of  calculations  based  on  these  equations  is 
presented. 
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II.    BACKGROOND  THEORY 

A.    PLASMA  PHYSICS 

L       Plasma    C ha racte ris tire 

A  plasma  is  formed  when  atoms  are  ionized  into  positive 
ions  and  negative  electrons  such  that  they  form  a  gas. 
However,  not  all  such  gases  are  considered  to  be  plasmas.  A 
more  complete  definition  of  a  plasma  has  been  given  by  Chen 
[Ref.    1]: 

"A     plasma     is     a     quaslneutral     gas     of    charged     and     neutral 

particles    which  exhibits   collective   behavior.  " 
"Collective    behavior"    means    that    the    behavior    of    a    small    region 
of   the    plasma    is    determined    not   only   by    conditions    in   the    plasma 
that     immediately     surrounds     it,     but    also     by     conditions     in     the 
plasma    far   away   from   it. 

A  plasma  can  be  characterized  by  two  quantities,  n,  the 
plasma  density  and  T,  the  plasma  temperature.  Since  the  plasma 
must  be  quasineutral  the  density  of  ions,  ni,  and  the  density  of 
electrons,  ne,  must  be  approximately  equal  Therefore,  the 
plasma   density   is    defined  as 

/I*  na  «  nt.  (  ec!n-  2-1  ) 

The  plasma  temperature  is  generally  given  in  units  of  energy 
which  corresponds  to  kT,  where  k  is  Boltzmann's  constant  and  T 
is     the     plasma      temperature     in     Kelvin.      For     a     gas      in     thermal 
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equilibrium,     which  follows     a    Maxwellian    distribution,     the    average 

energy       of      the  gas       particles       is       related      to       the       plasma 
temperature    by 

3  (  e<?n.  2.2  ) 

where    the    plasma    is    assumed   to    be    three    dimensional 
2.     Debye-  ShteJding 

When  a  plasma  is  subjected  to  an  external  electric  field 
the  electrons  and  ions  will  drift  in  opposite  directions.  As 
they  drift  apart,  they  produce  a  charge  separation  and  an 
internal  electric  field,  which  opposes  the  external  field.  The 
width  of  the  region,  over  which  the  charge  separation  must  occur 
to  balance  the  external  electric  field,  is  proportional  to  the 
thermal  energy  of  the  plasma  particles.  If  the  dimension  of  the 
bulk  plasma  is  greater  than  the  dimension  of  the  region  over 
which  this  charge  separation  occurs,  the  interior  of  the  plasma 
will  be    shielded   from   the   external  field. 

The  width  of  the  charge  separation  can  be  determined, 
approximately,  by  considering  a  plasma  where  the  ions  are  fixed 
in  space,  over  the  time  frame  of  interest,  and  the  electrons 
obey  a  Maxwellian  distribution  such  that  the  Boltzman  relation 
applies.      Using  Poisson's    Equation   in   one    dimension  gives 

d2V  e  (  A  C  eqn.  2.3  ) 
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where  V  is  the  electric  potential  and  e  is  the  electric  charge. 
Since  the  ions  are  considered  fixed,  the  density  of  ions  will  be 
constant  and  equal  to  the  plasma  density  no.  The  Boltzman 
relation  then  gives  the  electron  distribution  in  the  region  of 
interest   as 

f  eV  \  (  eqn.  2.4  ) 


n&  =  R0exp 


kT& 


Substituting  equation  2.4    into   equation  2.3    gives 

To  solve  equation  2.5  the  exponential  can  be  expanded  in  a 
Taylor  series,  and  neglecting  terms  of  second  order  and  higher 
gives 

d2V       ne2  C  eqn-  2.6  ) 

dx2     e0kTe 

The  solution  of  equation  2.6,  which  is  a  homogenous  second  order 
linear   ordinary   differential  equation   is    given   by 

(    \x\\  C  ec?n'  2-7  ) 

K  =  K0exp    -!—* 
\     ad  j 

where 


€0kT.  (eqn.  2.8) 

This    quantity    is    known   as    the    Debye    Length   and   is    a    measure    of 
the    thickness    of    the    shield    that    screens     the    bulk    plasma    from 
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the    affect  of   external  fields.      This    screening   applies    as    long   as 
the     dimensions     of     the     bulk    plasma     is     much     greater     than     the 
Debye    Length. 
3.      Sheaths 

When  a  plasma  comes  in  contact  with  a  wall,  the  electrons 
and  ions  that  hit  the  wall  will  recorabine.  But  since  the 
electrons  are  moving  at  a  higher  velocity  than  the  ions,  the 
electrons  will  be  lost  faster  than  the  ions  resulting  in  a  net 
positive  charge  at  the  plasma- wall  boundary.  The  wall  will  then 
be  at  a  lower  potential  than  the  bulk  plasma  and  an  electric 
field  will  exist.  Figure  2.1  illustrates  the  potential  variation  in 
a  plasma  which  is  in  contact  with  a  wall  Due  to  Debye  Shielding 
a  layer  of  charge  separation  will  exist,  called  the  sheath,  which 
will  isolate  the  bulk  plasma  from  this  potential  difference.  The 
effect  of  the  sheath  is  to  accelerate  the  ions  and  to 
decelerate  the  electrons  that  enter  the  region  until  the  flux  of 
ions  is  balanced  by  the  flux  of  electrons.  Therefore,  only  those 
electrons  in  the  high  energy  tail  of  the  velocity  distribution 
will  be  energetic  enough  to  cross  the  sheath  and  all  others  will 
be    repelled   back  into    the    plasma. 

The  variation  of  the  potential  in  the  plasma  sheath  will 
now  be  considered.  In  order  to  simplify  the  problem  the 
following  assumptions    will  be    made: 

*  The    ions    enter   the   sheath  region  with  a   drift  velocity 

Uo, 

*  The    ions    have    Ti=0    so    that  all  ions    have    a   velocity   uo 
at   the    plasma-sheath   boundary, 
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*  The   sheath  region  is    collisionless    and   in  steady   state, 

*  The    potential  decreases    monotonically   with  x, 

*  The   electrons    follow   the   Boltzman  relation  (equation 
2.4). 


0  =  0 


•©     © 


■»-  x 


Figure    2.1      Plasma-Wall  Interaction 


Applying   conservation   of   energy    to    the    ions    in   the    sheath 
region   gives 


l-Mu2-l-Mu20-eV(x). 


(  eqn.  2.9  ) 


Solving  equation  2.9    for   the    ion   velocity   gives 


a-  Ju.1- 


2eV 
Kf 


(  eqn.  2.10  ) 


The     ion    density    within    the    sheath    is     determined    from    the     ion 
equation  of  continuity   which  is    given   by 

at{x)iL{x)  -  n0u0.  (  eqn.  2.11   ) 
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Using    equation    2.10     for    the     ion    velocity,     equation    2.11    can    be 
solved   for   the    ion  density   which  gives 

1  (eqn.2.12) 


ni{x)  =  n0[  1  - 


2eV 


Substituting  equation   2.4    and   2.12    into   equation  2.3    gives 

(  eqn.  2.13  ) 


dzV     en0 


dx 


exp 


V\ 


kTj 


-     1- 


2eV 
u20M 


Equation  2.13  is  the  non-linear  planar  sheath  equation.  By 
multiplying  each  side  by  the  first  derivative  of  the  potential 
the    equation  can   be    integrated   once    to   give 

V'2-V'l  = 


2n0kT 


exp| 


qV 
kTe 


-  1 


kTa 


1- 


2eV 


-  1 


(  eqn.  2.14  ) 
Further  solution  of  equation  2.14  would  require  numerical 
methods.  If  the  electric  field  inside  the  plasma  is  zero,  then 
the  first  derivative  of  the  potential  inside  the  plasma  must 
also  be  zero.  The  left  hand  side  of  equation  2.14  is  therefore 
greater   than  zero,    and   the    following   inequality   results 

(  eqn.  2.15  ) 


(  eV  ]      ,      Mill 
expl  r^  I-  1  + 


kT9 


kT, 


1  - 


2eV 


-  1 


>0. 


16 


This  inequality  can  be  greatly  simplified  by  expanding  the  left 
hand  terms  in  Taylor  series  expansions  and  neglecting  terms  of 
third  order  and  higher.  The  resulting  inequality  is  known  as  the 
Bohm   sheath  criterion  and   is    given   by 


tin   > 


kT e  C  ec*n-  2-16  ) 


M 

This  requires  that  the  ions  must  enter  the  sheath  region  with  a 
velocity  greater  than  the  ion  acoustic  velocity  of  the  plasma. 
In  order  for  this  to  occur,  there  must  be  a  finite  electric  field 
within  the  plasma.  Therefore,  the  assumption  that  the  potential 
and  the  electric  field  are  zero  at  the  plasma-sheath  boundary  is 
only   an   approximation. 

B.        SPACE    CHARGE    LIMITED    CURRENT    FLOW 
L     CMMzLangnuir  Laja. 

The  current  that  can  flow  between  the  cathode  and  anode 
of  a  vacuum  diode  is  limited  by  the  space  charge  of  the 
electrons  that  exist  between  them.  The  electrons  distort  the 
external  field  and  thereby  reduce  and  can  even  reverse  the  field 
at   the    cathode    surface. 

As  an  example,  consider  a  diode  consisting  of  infinite  flat 
plates  separated  by  a  distance  d.  Assume  that  there  is  an 
unlimited  supply  of  electrons  available  from  the  cathode,  and 
that  the  initial  velocity  of  the  electrons  after  emission  is 
negligible    compared    with    the     velocity    that    they    will    gain    while 
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crossing  the  electrode  gap.  The  kinetic  energy  of  the  electrons 
crossing  the  gap  and  the  current  density  carried  by  the 
electrons    is    given   by 

l-mv2~eV  (  eqn.  2.17  ) 

and 

j  -  env .  (  eqn.  2.18  ) 

Using  Poisson's  equation,  in  one  dimension,  the  current- voltage 
relationship  can   be    derived   as    follows: 

d2V(x)     &n(x) 

,    z     =— •  (eqn.  2.19) 

dx*  e0 

The    electron  velocity    is    determined   from  equation  2.17    to    be 

v(x)-J  —  V\  (  eqn.  2.20  ) 

V  m 

The  electron  density  is  determined  using  equation  2.18  and  2.20, 
and   after   substituting  into   equation  2.19    gives 

p-2.J-Mv-K  (eqn.  2.21) 

dx*     e0  V  2e 

In  order  to  solve  equation  2.21  both  sides  are  multiplied  by  the 
first  derivative  of  the  potential  and  reduced  to  the  following 
form 

d  (dVY     j     l2m    -\dV 

V  '=— .  (  eqn.  2.22  ) 


(dVV m  j     12m 
\dx)      e0V    e 


dxKdxj       €0yj    e  dx 
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Equation  2.22  is  then  integrated  over  the  gap  using  the 
boundary  conditions  that  the  electric  field  and  the  potential  are 
zero    at   the    cathode.      The    result   of   this    integration   gives 


dV  2j    /2m,, 

—  -J— J V  .  (  eqn.  2.23  ) 

dx      V  eo  V    e 

Integrating  once  again  over  the  gap,  applying  the  boundary 
conditions,  and  rearranging  gives  the  desired  current- voltage 
relationship  of 

4/2  <f0     fjv} 

9       V  rnx 


Jm'-^—\l-Zz'  C  ec?n-  2-24  ) 


Equation  2.24  is  known  as  the  Child-Langmuir  Law  for 
space  charge  limited  current  flow  [Ref.  2].  Although  it  has  been 
derived  here  using  simple  assumptions,  the  proportionality  of  the 
current  density  to  the  3/2  power  of  the  potential  difference 
remains  for  more  difficult  geometries  under  non-relativistic 
conditions.  It  should  also  be  pointed  out  that  the  Child-Lang- 
muir Law  applies  to  both  positive  and  negative  charge  carriers. 
Figure  2.2  is  a  plot  of  electron  current  density  versus 
potential  for   several  gap  widths. 

2.     Space.  Charge  r.imit^H  Current  in  a  Sheath 

In  deriving  the  planar  sheath  equation  it  was  assumed 
that  the  electron  density  followed  the  Boltzman  relation 
equation  2.4.  Therefore,  in  a  region  close  to  the  wall  the 
electron  density   will  be    much   less    than  the    ion  density 
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Figure    2.2      Child- Langmuir   Space    Charge    Limited   Current   Flow 
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and  can  be  neglected.  The  first  derivative  of  the  potential  at 
the  point  in  the  sheath  where  the  electron  density  can  be 
neglected,  is  also  assumed  to  be  negligible  compared  with  the 
slope   of  the    potential  near   the    wall      Equation  2.13    then   reduces 

to 

d2V        en0(       2eV)"i 

— -i--— -    1--TT7       •  (  eqn.  2.25  ) 

Making   the    change    of   variable    that   V  =    -V   and   noting   that 

-T-  »  1  (  eqn.  2.26  ) 

equation  2.25    reduces    to 
dx2        e0     V  2e 

The  ion  current  density  is  constant  across  the  sheath  and  is 
given   by 

ji  =  en0u.0.  (  eqn.  2.28  ) 

After      substituting      equation      2.28  into       equation      2.27,       the 

resulting    expression    is     the     same     as  equation    2.21.      Therefore, 

the     solution     of     equation     2.27      for  the      ion     current     density 

results  in  the  Child-Langmuir  Law  for  space  charge  limited 
current   flow,    equation   2.24. 
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3.       MacKeocm    Equation 

In  an  arc,  the  current  is  carried  primarily  by  the 
electrons  emitted  by  the  cathode  and  by  a  smaller  amount  of 
positive  ions  flowing  from  the  plasma  to  the  cathode.  If  the 
space  charge,  which  causes  the  cathode  drop  in  the  plasma 
sheath,  is  high  enough  it  may  exert  a  strong  electric  field  at 
the  surface  of  the  cathode.  This  can  then  result  in  a  large 
field   emission  current  in  addition  to   any   thermionic   current. 

The  electric  field  at  the  cathode  surface  can  be 
determined  using  Poisson's  Equation  The  following  derivation 
measures  distance  and  potential  difference  from  inside  the 
plasma  where  the  cathode  drop  begins.  The  electron  and  ion 
current   densities    are    given   by 

jQ  =  QnaV  (  eqn.  2.29  ) 

and 

y*i  =  entu.  (  eqn.  2.30  ) 

Substituting  equation   2.29    and   2.30    into    equation   2.3    gives 

d^V=s_\_fj1_h\  C  eqn.  2.31   ) 

dx2        €0\  a      v  J 

where  V  is  the  potential  difference,  ji,  the  ion  current  density, 
je,  the  electron  current  density,  u,  the  ion  velocity,  and  v,  the 
electron  velocity. 

The  velocity  of  the  ions  and  electrons  are  determined 
using     the     conservation     of     energy.      The      initial     energies      are 
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neglected,  based  on  the  assumption  that  the  energy  gained  by 
the  ions  and  "electrons  in  crossing  the  sheath  are  much  greater 
than  the    initial  energies.      For   the    ions    the    energy    balance    is 


-Mu2  =  eV 
2 


and   for   the    electrons    the    energy    balance    gives 


(  eqn.  2.32  ) 


l-mv2  =  e{Vc-v) 


(  eqn.  2.33  ) 


where  M  is  the  mass  of  the  positive  ion,  ra,  the  mass  of  the 
electron  and  Vc  is  the  potential  of  the  cathode.  Substituting 
equation  2.32    and   2.33    into    equation  2.31   gives 

(  eqn.  2.34  ) 


d/V_ 
dx2 


1 


M 


\2eV 


J. 


m 


2e{Vc-V) 


Multiplying    both    sides     by    the    first    derivative    of    the    potential 
difference,    equation  2.30    can   be    integrated   once    and   gives 


'-WZ 


+  J.i 


m[Vc-V) 


m 


V,\ 


2e 


y* 


2e    J 


(  eqn.  2.35  ) 
where  F  is  the  electric  field.  Evaluating  for  V  -Vc  gives  the 
electric   field   at   the    cathode    surface    as 

(  eqn.  2.36  ) 
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Equation  2.36  is  known  as  Mackeown's  Equation  [Ref.  3]  and 
can  be  used  to  determine  the  electric  field  existing  at  the 
cathode  provided  the  values  of  the  cathode  drop,  the  ion 
current  density,  and  the  electron  current  density  are  known. 
Figure  2.3  is  a  plot  of  electric  field  at  the  cathode  surface 
versus  ion  current  density  for  several  ratios  of  electron  to 
ion  current   densities. 

C.    ELECTRON  EMISSION  FROM  METALS 

The  theory  of  electron  emission  from  metals  is  very  well 
understood.  The  basic  processes  have  been  explained  on  the 
basis  of  a  free  electron  gas  in  the  metal,  that  is  described  by 
Fermi- Dirac  statistics,  and  a  potential  barrier  at  the  metal 
surface.  Figure  2.4  shows  a  typical  potential  energy  diagram 
for  an  electron  as  a  function  of  distance  from  the  metal- vacuum 
interface.  For  thermionic  emission,  that  is  emission  at  high 
temperature  and  low  electric  field,  the  electrons  are  emitted 
predominately  by  passing  over  the  potential  barrier  at  the  metal 
surface.  For  field  emission,  that  is  emission  at  low  temperature 
and  high  electric  field,  the  electrons  are  emitted  predominately 
by   tunneling   through  the    potential  barrier. 
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ELECTRIC  FIELD  AT  CATHODE  SURFACE 
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Figure    2.3      MacKeown's    Equation   for   Electric   Field  at   the 

Cathode   Surface 
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VACUUM 


Figure    2.4      Potential  Energy   Diagram   for  a   Metal- Vacuum 

Interface 

In  the  past,  thermionic  and  field  emission  have  been  treated 
separately,  which  has  resulted  in  the  well  known  Richardson 
equation  for  thermionic  emission,  and  the  Fowler-  Nordheim 
equation  for  field  emission.  The  expressions  that  will  be  used 
here,  however,  are  based  on  an  approach  that  treats  both 
thermionic  and  field  emission  from  a  unified  point  of  view 
[Ref.    4]. 

Hartree  units  are  used  in  the  following  sections  because 
they  greatly  simplify  the  expressions.  Thus  the  following 
quantities    are    redefined   as: 


J  -; 


m2e9 


{4neoyK7, 


(  eqn.  2.37  ) 
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m2e' 


F-f(^OT'  (eqn.  2.38) 


me* 


(fcT)'»(fcT)- — ,  (  eqn.  2.39  ) 

(4/re0r£2 

m  s 

t'**r. ia«'  (eqn.  2.40) 

(4/re0)2fi:2 

where  the  starred  quantities  on  the  left  are  in  SI  units  and 
their  respective  counterparts  on  the  right  are  in  Hartree  units. 
The  following  sections  will  give  expressions  that  can  be  used 
to  determine  the  current  density  that  can  be  expected  to  be 
emitted  from  a  metal  given  the  surface  temperature,  the  electric 
field  at  the  surface,  and  the  work  function  for  the  metal  under 
consideration.  In  addition,  a  set  of  applicability  conditions  are 
given  that  can  be  used  to  determine  the  validity  of  the 
particular  expression  for  a  given  situation.  Figure  2.5 
illustrates  the  regions  of  applicability  of  the  emission 
equations  on  a  plot  of  temperature  versus  electric  field. 
Appendix  A  is  a  Fortran  program  that  calculates  the  current 
density  that  can  be  expected  to  be  emitted.  Appendix  B  and  C 
are  Fortran  subroutines  that  calculate  the  current  density  and 
applicability   of  the   expressions    to    be    used   respectively. 

L        T>iftT-mir>nir!-gi^1H    KmiftSJOT) 

The    current    density    that    can    be    expected    from    thermionic 
emission,    corrected   for   the    Schottky   effect,    is    given   by 
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(  eqn.  2.41   ) 
(  eqn.  2.42  ) 
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Figure   2.5      Regions    of   Applicability   for   Thermionic,    Field,    and 
Intermediate    Electron  Emission   Expressions 

For    equation    2.41    to    be    applicable    the    following    conditions    must 
be    satisfied*. 

(l-d)\  1  _-%(  _\\  (eqn.  2.43) 


In 


d(l  -d 


>-jiF  4   0-F2    , 
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/(l-d)\      1  -I 

lnr   ..-■       )-- ->-jzF   \ 


(  eqn.  2.44  ) 


d      j     [l~d) 

Equation   2.41   reduces    to    the    Richardson   equation   in   the    limit    as 

the     electric     field     is     reduced     toward     zero.      Appendix     D  is     a 
listing        of        a        Fortran        subroutine        that        calculates        the 

Thermionic- Field   emission  current   density    based   on   equation  2.41. 
2.        Field- Th^TTninnin    EIniSfiiQI] 

The      current     density      that     can     be      expected     from  field 
emission,    corrected   for   temperature,    is    given   by 


J  = 


rtckT 


1  6/r20t2\sin  TtckT 


exp     - 


• 

4/2Vu 


where 


t(y)  = 


1+y 


(l+y)£ 


1  -y 
1+y 


3F 


-yK 


i-y 
l+y 


(  eqn.  2.45  ) 


(  eqn.  2.46  ) 


c  =  2/2^-t(y)( 


v(y)  =  -JJ-i-2E 


y-i 
2y 


+  (y+I)tf 


for     y  >  1  , 


Ky)=vT+y \E 


for     y  <  1  , 


1 

y  -  -r 


i-y 

I+y 


-yK 


y-i 
2y 


i-y 
l+y 


(  eqn.  2.47  ) 


(  eqn.  2.48  ) 


(  eqn.  2.49  ) 
(  eqn.  2.50  ) 
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Equations  2.46,  2.48,  and  2.49  are  the  Fowler- Nordheim  Elliptic 
functions.  The  functions  K[k]  and  E[k]  are  complete  elliptic 
integrals     of     the     first     and     second     kind     respectively     and     are 

defined   as 

« 
rl  -i 

K[k]~  I    (l  -fc2sin20)  3dOt  (eqn.2.51) 

Jo 

rl  I 

E[k]-       (l  -k2s\n29YdO.  (  eqn.  2.52  ) 

Jo 

For    equation    2.45    to    be    applicable    the    following    conditions    must 

be    satisfied: 


F  kT 

tf-Fi>—  +  r—T,  (  eqn.  2.53  ) 

Tt  [  -  CK  1 

\-ckT>i2]kTt  (  eqn.  2.54  ) 


J2<t>3 
/- — tt- ^y(y)-  (  eqn.  2.55  ) 

Equation  2.45  reduces  to  the  Fowler- Nordheim  equation  in  the 
limit  as  the  temperature  approaches  zero.  Appendix  E  is  a 
listing  of  a  Fortran  subroutine  that  calculates  the 
Field- Thermionic  current  density  based  on  equation  2.45.  This 
subroutine  makes  calls  to  the  subroutines  listed  in  Appendix  F, 
G,  and  H,  which  are  used  to  calculate  the  Complete  Elliptic 
Integral  of  the  First  Kind,  Complete  Elliptic  Integral  of  the 
Second  Kind,  and  the  Fowler- Nordheim  Elliptic  functions 
respectively. 
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3.        TnteiTDftriiata    Fift1ri-Thi»nirinTvic»    Rmissinn 

When  equations  2.41  and  2.45  are  not  applicable,  the 
current  density  from  fie  Id- thermionic  emission  may  be  determined 
in   the    intermediate    region   as 


F  fkTtV         f      <fi  F2G      \ 

J'^Zl  ^T=      exP + 


2n\2n  J       r\    kT     24(fcT)3; 
where 


(  eqn.  2.56  ) 


(9  =  — -.  (  eqn.  2.57  ) 

t2     t3 

The    argument   of   the    functions    t,    equation   2.46,    and    v,    equations 
2.48    and   2.49,    as    used   in   equations    2.56    and   2.57,    is    given   by 

y„~— —  .  (  eqn.  2.58  ) 

Eta     can     be     determined     by     solving     the     following     equation     by 
iteration: 

F2 
r?  = ——. r.  (  eqn.  2.59  ) 

8(fcT)2C2(yJ 
For    equation   2.56    to    be    applicable    the    following    conditions    must 
be    satisfied: 

->  1  +     ,.    °     .  (eqn.  2.60) 

2j2t[yn) 
dm . — — ,  (  eqn.  2.61  ) 


F 2 


>-<p  + — .  (  eqn.  2.62  ) 


8(fcT)2t"  l-F{2^kTtty' 
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Appendix     I     is     a     listing     of     a     Fortran     subroutine      that 
calculates      the      intermediate      Fie  Id- Thermionic     emission     current 
density   based   on  equation  2.56. 
4.     Explosive  Electron  Emission 

The  emission  mechanisms  discussed  so  far  have  wide 
applicability  to  many  physical  devices.  However,  with  the  advent 
of  high  current  cathodes  the  thermionic  and  field  emission 
theories  were  unable  to  explain  the  high  current  densities 
observed.  One  explanation  for  this  high  current  density  is  that 
the  surface  of  the  cathode  is  not  smooth,  but  actually  has  a 
surface  covered  with  micro  points,  or  whiskers.  The  significance 
of  the  whisker  is  that  the  electric  field,  at  the  surface  of  the 
whisker  is  enhanced  by  a  factor  of  10  -  1000.  Figure  2.6 
illustrates  the  electric  field  enhancement  that  can  occur  at  the 
tip  of  a  whisker.  This  results  in  a  higher  field  emission 
current  density  from  the  whisker,  and  ultimately  leads  to 
explosive  destruction  of  the  emission  site.  This  destruction 
produces  a  local  burst  of  plasma  that  expands  rapidly  away  from 
the  cathode  surface  and  has  therefore  come  to  be  known  as  a 
cathode  flare  or  plasma  jet.  The  electrons,  that  reach  the 
anode,  are  then  emitted  from  the  surface  of  the  plasma  jet. 
This  whole  process  is  commonly  referred  to  as  explosive  electron 
emission  [Refs.    5    and   6]. 
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Figure   2.6      Electric   Field   Enhancement   at   a    Whisker 


Fursei  and  Zhukov  [Fefs.  7,  8,  and  9]  have  experimentally 
characterized  the  explosive  emission  process  and  have  confirmed 
the  existence  of  raicropoints  and  that  the  current  density 
follows  the  Fowler-  Nordheim  field  emission  theory  with  the 
inclusion  of   a    field   enhancement   factor. 
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III.        UNTPOLAR    tRQJMQ. 

A.        Robson-Thonemann    Model 

In  1958,  Robson  and  Thoneman  observed  an  electrical  arc 
burning  between  a  dense  hot  plasma  and  a  metal  wall  [Ref.  10]. 
What  made  this  unusual  was  that  there  was  no  anode  present  to 
collect  the  electrons  emitted  from  the  conducting  wall 
Apparently,  the  wall  was  acting  as  both  a  cathode  and  anode,  and 
therefore,    they   named   what   they   observed  a   "unipolar   arc." 

Robson  and  Thoneman  proposed  the  following  model  to  explain 
the  occurrence  of  unipolar  arcing.  When  a  hot  dense  plasma 
comes  in  contact  with  a  conducting  wall  a  sheath  will  form.  If 
the  plasma  temperature  is  high,  the  sheath  or  floating  potential 
can  exceed  that  necessary  to  initiate  an  electric  arc.  A 
cathode  spot,  which  is  a  local  center  of  emission,  could  then 
form  on  the  metal  surface.  The  large  electron  current  flow  into 
the  plasma  from  the  arc  would  result  in  lowering  the  plasma 
potential  to  the  cathode  fall  potential  of  the  arc.  This  allows 
more  electrons  from  the  high  energy  tail  of  the  velocity 
distribution  to  cross  the  sheath  and  recombine  at  the  wall  and 
thus  closing  the  current  loop.  Figure  3.1  illustrates  the 
situation  in  the  sheath  prior  to  initiation  of  an  arc,  and  Figure 
3.2    shows    the    return  current   flow   after  arc   initiation. 


34 


Figure    3.1     Sheath  Region  Prior   to    Unipolar   Arc   Initiation 
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Figure    3.2      Sheath  Region   After   Unipolar   Arc   Initiation 


The    resulting    change    in    current    flow    outside    of    the    arc    is 
then   given   by 

(  eqn  .3.1    ) 


1=  Aen 


2nm 


exp 


exp| 


qV  f 
kf] 
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where  I  is  the  current,  A  is  the  area  of  the  metal  plate,  and  Vc 
is  the  cathode  drop  of  the  arc.  But  since  the  net  current  flow 
between  the  plasma  and  the  wall  must  be  zero,  equation  3.1  also 
gives    the   current  flow   inside   the   arc. 

The    ion    current    density    in    the    sheath    can    be    approximated 
as 


kj  (  eqn.  3.2  ) 

where  the  Bohm  sheath  criterion  is  used  to  approximate  the  ion 
velocity.  The  electron  current  density  in  the  sheath  can  be 
determined  by  noting  that  the  net  flux  of  electrons,  *,  from  the 
plasma    into    the   sheath  is    given  by 

1  (  eqn.  3.3  ) 

4>  =  -nQv. 
4 

Using  the  Boltzman  relation  for  the  electron  density  and  the 
expression  for  the  average  speed  of  particles  following  a 
Maxwellian   distribution  gives 


1  foVA        SkT,,  (eqn.  3.4) 

/'   =-e/inexp    L  / . 

J*     4      °         \kTjyj    nm 

The      net      current     density      across      the  sheath     must     be      zero. 

Therefore,      setting     equation     3.2      equal  to     equation     3.4,      and 
solving  for   the   sheath  potential  gives 

kT0       f     M     \  (eqn.  3.5) 

f      2e       \2nmJ 
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Robson  and  Thoneman  point  out  that  there  is  a  minimum 
current  necessary  to  maintain  an  arc  burning.  Therefore, 
according  to  equation  3.1,  if  the  plate  area  is  small  enough  a 
unipolar  arc  can  not  occur.  For  example,  using  copper  with  a 
plasma  temperature  of  4  eV,  a  cathode  fall  potential  of  15  V,  a 
plasma  density  of  1013  cm-3,  and  an  arc  current  of  100  amperes 
would  require  a  plate  area  of  115  cm2.  For  a  circular  area,  the 
radius    required   is    6.0  5    cm. 

B.   Schwirzke-Taylor  Model 

Schwirzke  and  Taylor  have  expanded  on  the  basic  concepts  of 
the  Robson-Thoneman  model  [Ref.  11].  They  propose  that  for  an 
arc  to  occur,  the  ion  density  above  the  cathode  spot  must 
increase  to  allow  for  the  arc  current  to  flow.  This  is  contrary 
to  the  Robson-Thoneman  model  which  assumes  a  constant  plasma 
density.  Once  a  cathode  spot  forms  the  desorption  of  neutral 
atoms  and  evaporation  of  the  metal  atoms  from  the  spot  is 
greatly  increased.  If  only  a  small  fraction  of  these  neutral 
atoms  become  ionized  in  the  sheath  the  plasma  density  will 
increase  above  the  cathode  spot.  Thus,  according  to  equation 
2.8  the  Debye  length  will  decrease  resulting  in  a  decrease  of 
the  sheath  width  above  the  cathode  spot.  The  electric  field 
above    the    spot    increases,     where    the    field    can    be    approximated 


as 


r    I       V  f        I T^,     f      M      \  (eqn.3.6) 

KD  \  znm  J 
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Schwirzke  [Refs.  12  and  13]  elaborates  further  on  the 
electric  fields  set  up  inside  the  high  density  plasma  which  forms 
above  the  spot.  The  fluid  equation  of  motion  of  the  electrons 
in  the    plasma    is    given  by 

n.m—  =  -erc8(f  +  i,8X£)-V/\+  —  j. 

The  time  derivative  can  be  neglected  over  the  time  span  of  an 
arc,  and  further  assuming  no  magnetic  field,  equation  3.7  can  be 
solved   for   the    electric   field,    which  gives 

_      7         1  C  eqn.  3.8  ) 

F-- V/V 

a    ena 

The  first  term  on  the  right  hand  side  of  equation  3.8  can  be 
neglected  in  a  weakly  ionized  plasma.  This  implies  that  the  mean 
free  path  length  for  collision  between  the  electrons  and  the 
neutral  atoms  is  much  greater  than  the  dimensions  of  the  region 
of  increased  density  above  the  cathode  spot.  Thus,  the 
increased  plasma  pressure  above  the  spot  creates  a  radial 
electric   field  given   by 

_  1  kT9  (eqn.3.9) 

F V/5,' Vna. 

The  affect  of  this  radial  electric  field  is  to  lower  the  plasma 
potential  in  a  ring  around  the  cathode  spot.  The  change  in  the 
potential  can  be  determined  by  integrating  equation  3.9  from  the 
center  of  the  arc,  where  the  plasma  density  is  highest,  radially 
outward    until    the    plasma    density    equals    the    background    density. 
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The    change    in  potential  is    then  given   by 


o 


JK(r)- -In 

e       u(0y 


(  eqn.  3.10  ) 


In  this  ring  of  reduced  potential,  electrons  in  the  tail  of  the 
Maxwellian  distribution  can  cross  the  sheath  more  easily  and 
thus  a  higher  return  current  flow  from  the  plasma  to  the  wall 
results.  This  current  flow  closes  the  current  loop  of  the 
unipolar  arc.  The  concepts  of  the  Schwirzke- Taylor  model  are 
illustrated   in   Figure    3.3. 
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Figure    3.3      Schwirzke- Taylor   Unipolar   Arc    Model 
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If  equation  3.10  is  set  equal  to  the  floating  sheath 
potential,  equation  3.5,  the  ratio  of  the  electron  density  in  the 
plasma  above  the  spot  to  the  background  plasma  density  which 
reduces    the   sheath  potential  to   zero   can  be   determined   as 


n  M  (  eqn.  3.1  1    ) 


na0      v  2/rm 

Equation  3.3  and  equation  3.11  can  be  used  to  calculate  the 
maximum  area  necessary  to  support  a  given  unipolar  arc  return 
current.  By  setting  the  current  equal  to  the  current  derived 
for  the  Robson-Thoneman  model,  equation  3.1,  the  ratio  of  the 
return  current  areas  required  by  each  model  can  be  compared. 
The    result  of  the    comparison  gives 


AR.T    rari     r  *vc\    [j^v1    rw      (eqn. 3.12) 

^_r~V2/rm\eXPV"it7j"V     M    J     >  V  2~^ 

where  Ar-t  and  As  -  T  are  the  areas  required  by  the 
Robson-Thoneman  model  and  Schwirzke- Taylor  model  respectively. 
As  an  example,  for  copper  with  a  cathode  fall  potential  of  15 
volts  and  a  plasma  temperature  of  4  eV  the  area  ratio  is  8.43  x 
103.  Therefore,  for  a  100  ampere  arc  the  area  required  by  the 
Schwirzke- Taylor  model  is  only  1.37  x  10-2  cm2,  or  for  a  circular 
area    the    radius    required   is    only   6.59    x  10" 2    cm. 
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IV.   CATHODE  &ZQXS- 

A.    TYPES  OF  CATHODE  SPOTS 

A  cathode  spot  is  a  small,  luminous  region  on  the  cathode, 
where  a  localized  transfer  of  current  occurs  between  the 
cathode  and  a  plasma  in  an  arc  discharge.  The  type  of  arc 
discharge  can  be  either  bipolar  or  unipolar,  however  in  all  cases 
a  cathode  spot  is  present.  In  addition,  the  spots  are 
associated  with  the  erosion  of  the  cathode  and  therefore,  are 
also  localized  centers  for  the  transfer  of  material  from  the 
cathode  into  the  plasma.  This  has  been  evidenced  by  the 
formation  of   craters    on   the    surface    of  the    cathode. 

It  has  been  possible  over  many  observations  to  distinguish 
between  two  different  types  of  cathode  spots  [Ref.  14].  The 
first  type,  or  type  I  spots,  are  characterized  as  rapidly  moving 
spots  with  speeds  of  approximately  10 4  cm/s  across  the  surface 
of  the  cathode  and  with  relatively  low  erosion  from  the  surface. 
The  type  II  spots,  on  the  other  hand,  have  speeds  of 
approximately  10 2  cm/s  and  have  much  greater  erosion  from  the 
surface.  The  type  EC  spot  also  has  a  tendency  to  group  together 
and  form  clusters.  The  mechanism  of  erosion  for  the  type  I 
spot  is  believed  to  be  due  to  no  n- thermal  means,  such  as 
explosions  of  whiskers  or  micro  in  ho  mo  genie  ties  on  the  cathode 
surface.      The    erosion  mechanism   for   the    type    II   spot    is    believed 
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to    be    due    to    thermal  evaporation. 

B.    THKOHY  OF  STATIONARY  CATHODE  SPOTS 

The  theoretical  modeling  of  cathode  spots  has  progressed 
over  the  last  100  years,  but  up  until  the  present  time  there 
has  been  no  satisfactory  understanding  of  the  processes 
occurring  within  the  spot.  The  difficulty  occurs  due  to  the 
diverse  phenomena,  which  in  themselves  are  extremely  complicated 
to  describe  mathematically,  that  take  place  within  the  small 
cathode- plasma  region  known  as  a  cathode  spot.  A  detailed 
analysis  would  require  an  extensive  use  of  the  kinetic  theories 
of  Plasma  Physics,  Solid  State  Physics,  Physics  of  High 
Temperature  Phenomena,  and  Gas  Dynamics.  In  order  to  avoid 
some  of  the  complications  the  majority  of  models  developed  for 
the  cathode  spot  are  stationary,  and  therefore  are  unable  to 
explain  the  initiation,  motion,  and  disappearance  of  the  cathode 
spot  [Refs.    12    and   15]. 

The  energy  balance  method  developed  by  Lee  and  Greenwood 
[Ref.  16]  has  been  the  basis  of  the  majority  of  stationary 
models  including  those  proposed  by  Kulyapin  [Ref.  17],  Kubono 
[Ref.  18],  Beilis  [Ref.  19],  and  Moizhes  and  Nemchinskii  [Refs.  20, 
21,  22,  and  23].  The  latter  work  by  Moizhes  and  Nemchinskii  are 
the  basis  for  the  model  developed  in  this  paper.  The  cathode 
spot   is    divided   into   three    regions    as    illustrated   in   Figure    4.1. 
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Figure    4.1      Model  of   the    Cathode-Plasma    Region 

Region  I,  referred  to  as  the  Langmuir  region,  is  a  thin 
space  charge  sheath  over  which  the  cathode  potential  drop,  Vk, 
occurs.  It  is  in  this  region  that  the  electrons  emitted  by  the 
cathode  are  accelerated  into  the  plasma  and  where  the  ions 
produced  in  the  generation  region  are  accelerated  towards  the 
cathode.  Region  II,  referred  to  as  the  generation  region  or 
Knudsen  Region,  is  where  the  neutral  atoms  evaporated  from  the 
cathode  are  ionized.  Region  III  is  an  expanding  channel  or 
transition  region  over  which  the  highly  dense  current  channel  in 
regions    I  and   II  expands    out   into    the    bulk  plasma. 
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In  developing  a  theoretical  model  of  a  cathode  spot,  one 
of  the  first  decisions  that  must  be  made  is  the  electron 
emission  mechanism.  In  general,  most  models  use  a  combination 
of   thermal  and   field   equations    that  are    expressed   as 

j9  =  f{Tk,F)  (eqn.4.1   ) 

where     je     is     the     electron    current    density,     Tk     is     the     cathode 

temperature,    and    F    is    the    electric    field   at   the    cathode    surface. 

In    actual    use,     expressions    such    as    equation    2.41,     2.45,     or    2.56 

are       used.      The       particular       equation      used       depends       on       the 

applicability    under    the     given    conditions.      The    electric    field    at 

the     cathode     surface     is     determined     using     MacKeown's     Equation, 

equation   2.36. 

An    energy    balance    must    exist    in    the    plasma    region    above 

the    cathode    spot   and   can   be    described   by   equation   4.2. 


^{eVk  +  2kTk)-J^{Et  +  2kT,)  +  fi^(ja  +  jl) 


*.T_DI  .        .  ,  (  eqn.  4.2  ) 

e 


Energy  is  brought  into  the  region  by  the  electrons  emitted  from 
the  cathode  and  consists  of  the  thermal  energy  and  the  energy 
gained  as  the  electrons  were  accelerated  through  the  cathode 
potential  drop.  Energy  is  removed  from  the  region  by  ions 
flowing  back  to  the  cathode.  This  energy  loss  consists  of  the 
ionization  energy  and  the  thermal  energy  of  the  ions.  A  further 
loss  occurs  from  both  ions  and  electrons  that  leave  the  region 
and    pass     into     the    bulk    plasma    and    is    described    by    the     Peltier 
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heat  flux  ;0kT,.  In  a  plasma,  in  which  Coulomb  collisions  are  the 
predominant  scattering  process,  the  energy  transfer  coefficient, 
0,  is  approximately  equal  to  3.21.  It  should  be  noted  that 
electrons  in  the  plasma,  also,  can  leave  the  region  towards  the 
cathode,  but  only  if  energetic  enough  to  over  come  the  cathode 
potential  drop,  Vk.  This  reverse  electron  current  can  be 
approximated   by 


erjaP 


f     QV„\  (  eqn.  4.3  ) 


rev 


exp 


For  low  plasma  temperatures  and  high  cathode  potential  drops 
the  reverse  electron  current  density  is  much  less  than  the  ion 
current   density   and   can   therefore    be    neglected   in  this    model 

An   energy    balance    at    the    cathode    surface    must   also    exist 
and   can   be    described   by   equation   4.4. 


4KRTk  =  JtR2 


!+{Ei  +  2kTp  +  9Vk-t)-!j{+  +  2kTk)-gEvap 


(  eqn .  4.4  ) 

The  predominant  loss  mechanisms  considered  are  heat  conduction, 
emissive  cooling,  and  evaporation.  The  primary  energy  source 
considered  is  ion  bombardment  of  the  cathode  surface.  The  ions 
bring  to  the  surface  thermal  energy,  ionization  energy,  and 
energy  gained  while  being  accelerated  by  the  cathode  potential 
drop.  Energy  is  released,  approximately  equal  to  the  work 
function,  when  ions  and  electrons  recombine  at  the  surface.  The 
electrons    emitted    by    the    surface    remove    thermal   energy    and    an 
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energy  approximately  equal  to  the  work  function.  In  addition, 
neutral  atoms,  evaporating  from  the  surface,  remove  energy 
equivalent  to  the  heat  of  vaporization.  Heat  conduction  from 
the  cathode  spot  is  based  on  a  circular  spot  of  radius  R  at  a 
constant  temperature  Tk  and  where  k  is  the  thermal  conductivity 
of  the  cathode.  The  work  function  to  be  used  in  equation  4.4  is 
corrected   for   the   Schottky   effect  and   is    given  by 


&3F 

4/re0 


0  =  00- a/  tztt-  (  ecin'  4-5  ) 


The    total  current   flow   in  the    arc   is    given   by 

/  =  7tR2[ja+  jt).  (  6<*n-  4-6  ) 

A    fraction    of    the     neutral    atoms     evaporated    from    the     cathode 

will     return     to      the     cathode.      The     return     coefficient     can     be 

determined    from    equation    4.7    which    is     based    on    a    Monte    Carlo 

calculation   performed   by   Nerachinskii  [Ref.    24]. 


The    ion  current  density   is    determined  from 

/i«  -enJL  (  eqn.  4.8  ) 

4 


where 


nt-— —  (  eqn.  4.9  ) 


and 
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a 


8kTk 


(  eqn.  4.10  ) 


Substituting    equations     4.9     and    4.10     into    equation    4.8    gives     the 
ion  current   density   as 
eP 


Jt  = 


■\2nUkT 


n*{T,) 


(  eqn .  4. 1  1   ) 


where  P  is  the  vapor  pressure  of  the  cathode  material  and  «  is 
the  degree  of  ionization  The  degree  of  ionization  is  determined 
from   Saha's    Equation   which  is    expressed   as 


a  = 


3 


kTp\M kT  p 


P    (  2nh2  Y         f   Et 
exp 


kT 


(  eqn.  4.12  ) 


The    erosion   rate    of   neutral  atoms    is    given   by 
P 


g  = 


■\2nUk7 \ 


O-r?) 


(  eqn .  4. 1 3  ) 


The      potential      drop      inside      the      quasineutral      plasma      of      the 
transition   region   is    given   by 

/ 


U- 


\2noR)\  j 


^](l+z)-/?KT, 


ly] 


(  eqn.  4.1  4  ) 


where    z    is    the    average    charge    of    the    ions    in    the    region,     j    is 
the    total  current   density   given   by 


and   o   is    the    electric   conductivity   given   by 


(  eqn .  4.15  ) 
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1  .53x  10'2     '- 
a — T\  (eqn.4.16) 

where    the   Coulomb  Logarithm   is    approximated  as 


(    1       {P\ 

nA  =  23-ln 


(  eqn.  4.17  ) 


\TiMkJ 

The    total  potential  drop  across    the   arc   is    given   by 

V*rc  =  Vk  +  U-  (  eqn.  4.18  ) 

2.      Method   Q±  Solution 

The  system  of  equations  presented  in  the  previous  section 
is  not  closed.  Therefore,  more  information  is  required  in  order 
to  solve  the  system.  This  will  be  accomplished  here  using  the 
Steenbeck  Minimum  Principle  [Refs.  25  and  26].  Simply  stated  this 
requires  that  the  total  potential  drop  across  the  arc  must  be  a 
minimum  for  a  stationary  arc  to  exist.  Given  an  arc  current  and 
cathode  temperature  the  system  of  equations  is  solved  for  the 
cathode  potential  drop,  plasma  temperature,  ion  current  density, 
electron  current  density,  erosion  rate,  potential  drop  in  the 
transition  region,  and  total  arc  potential  drop.  The  cathode 
temperature  that  defines  the  stationary  cathode  spot,  for  the 
given  arc  current,  is  that  cathode  temperature  that  minimizes 
the  total  potential  drop.  If  a  minimum  total  arc  potential  can 
not  be  found  for  a  given  arc  current,  then  a  stable  arc  can  not 
exist   for  this    current. 
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The  first  step  in  the  solution  is  to  determine  the  plasma 
temperature,  Tp,  and  cathode  drop,  Vk.  To  accomplish  this  a 
guess  is  made  for  the  plasma  temperature.  The  partial  pressure 
of   the    cathode    material  is    determined   using  [Ref.    27] 

log/5-—  +B\qqTp  +  CTp  +  D  (eqn.4.19) 

where  for  copper  A=  -1.752x104,  B=  -1.21,  C=  0,  and  D=  13.21.  The 
degree  of  ionization  is  determined  using  equation  4.12  and  the 
return  coefficient  is  determined  using  equation  4.6.  The  ion 
current  density  can  now  be  determined  using  equation  4.1L 
Experimentally,  it  has  been  observed  that  the  cathode  drop  is 
approximately  15  volts,  so  this  will  be  used  as  a  first 
approximation  in  order  to  determine  the  electric  field  at  the 
cathode  surface  using  equation  2.36.  The  Schottky  corrected 
work  function  is  determined  using  equation  4.5  and  the  electron 
current  density  is  determined  using  equation  4.1.  The  spot 
radius  can  now  be  determined  using  equation  4.6.  The  erosion 
rate  is  determined  from  equation  4.13.  Equation  4.2  and  4.4  are 
each  solved   separately   for   the    cathode    potential  which  gives 


i.M-i{£(£l+2*r.)  +  ,*r.(i  +  £) 


-2*7* 


(  eqn.  4.20  ) 


and 
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vk2'-^-T^-(J-U+2kTlc)-Ei-2kTD^<l>)  +  ^EvaD. 
k2     nRji     k     e\;'ivr  *)        t  P     rj     ^     vaP 

(  eqn.  4.21  ) 
The  potential  is  then  determined  using  both  equation  4.20  and 
4.21.  The  results  are  compared  and  the  plasma  temperature  is 
then  adjusted.  This  procedure  is  then  repeated  until  the 
potentials  determined  from  both  equations  are  in  agreement.  The 
conductivity  is  determined  using  equation  4.16  and  4.17.  The 
potential  drop  in  the  transition  region  is  then  determined  using 
equation  4.14.  Finally,  the  total  arc  potential  is  determined 
from  equation  4.18.  Appendix  J  is  a  Fortran  program,  which  uses 
the  method  described  above,  to  calculate  the  cathode  spot 
parameters. 
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The  Fortran  program  listed  in  Appendix  J  was  used  to 
calculate  the  parameters  of  the  cathode  spot  for  arc  currents 
of  100,  150,  175,  200,  300  and  400  Amperes.  The  program 
required  knowledge  of  several  material  parameters  which  are 
specified   in  Table    1. 


TABLE     1 
CATHODE    SPOT    MODEL    CONSTANTS 

Constant 

Value 

Material 

Copper 

Ion   Mass    (gr/mole) 

63.54 

Heat   of   Vaporization   (kilo joules) 

300.3 

Ionization   Energy   (eV) 

7.726 

Work  Function  (eV) 

4.4 

Thermal  Conductivity   (W/m-K) 

40 

Figure  5.1  illustrates  the  variation  of  the  total  arc 
potential  with  cathode  temperature  for  different  arc  currents. 
The  cathode  temperature  for  a  stable  cathode  spot,  and  thus 
the  parameters  of  the  cathode  spot,  is  determined  from  the 
minimum  of  the  total  arc  potential  It  can  be  seen  from  Figure 
5.1  that  a  stable  arc  does  not  occur  for  arc  currents  below  200 
Amperes. 
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Figure   5.1      Cathode    Spot   Stability 
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The  results  of  the  calculations  are  listed  in  Table  2  through 
7  for  arc  currents  of  100,  150,  175,  200,  300,  and  400  Amperes 
respectively.  In  addition  Figures  5.2  through  5.7  illustrate  the 
variation  of  the  cathode  potential  transition  region  potential 
and   total  arc    potential  for   each  arc   current. 

The  cathode  spot  parameters  for  stable  arc  currents  of  200, 
300,  and  400  Amperes  are  summarized  in  Table  8.  The  results 
are  in  general  order  of  magnitude,  agreement  with  experimental 
data  [Ref.  23],  however  stability  has  been  indicated  at  arc 
currents  between  100  and  200  Amperes.  This  discrepancy  occurs 
since  the  material  constants  listed  in  Table  1  are  not  really 
constants  at  these  temperatures.  A  more  detailed  comparison  is 
not  warranted,  since  the  published  experimental  data  for  the 
cathode  spot  parameters  vary  widely.  Confirmation  of  the  models 
of  the  cathode  phenomena  is  greatly  limited  by  this  lack  of 
consistent   experimental  results. 
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TABLE    2 

CATHODE    SP( 

100    AMPERE    ARC    ON    A 

COPPER    CATHODE 

Cathode 

Plasma 

Cathode 

Transition 

Total    Arc 

Temperature 

Temperature 

Potential 

Potential 

Potential 

(Kelvin) 

(eV) 

(Volts) 

(Volts) 

(Volts) 

3500 

1.49 

22.89 

8.72 

31.61 

3550 

1.54 

21.54 

9.21 

30.75 

3600 

1.59 

20.34 

9.70 

30.04 

3650 

1.65 

19.27 

10.18 

29.45 

3700 

1.70 

18.32 

10.65 

28.97 

3750 

1.76 

17.49 

11.10 

28.59 

3800 

1.82 

16.75 

11.53 

28.28 

3850 

1.89 

16.10 

11.92 

28.02 

3900 

1.96 

15.54 

12.26 

27.80 

3950 

2.03 

15.06 

12.53 

27.59 

4000 

2.11 

14.66 

12.71 

27.37 

4050 

2.20 

14.34 

12.78 

27.11 

4100 

2.29 

14.09 

12.70 

26.79 

4150 

2.39 

13.91 

12.44 

26.35 

4200 

2.51 

13.82 

11.94 

25.76 

Cathode 

Partial 

Degree    o-f 

Return 

Spot 

Temperature 

Pressure 

Ionization 

Coe-f-ficient 

Radius 

(Kelvin) 

(atm.) 

(10-*    m) 

3500 

10.9 

0.1625 

0.4726 

4.807 

3550 

12.6 

0.1712 

0.4751 

4.258 

3600 

14.5 

0.1803 

0.4777 

3.777 

3650 

16.6 

0.1897 

0.4804 

3.354 

3700 

19.0 

0.1996 

0.4831 

2.981 

3750 

21.6 

0.2101 

0.4859 

2.652 

3800 

24.5 

0.2213 

0.4888 

2.360 

3850 

27.6 

0.2332 

0.4918 

2.100 

3900 

31.1 

0.2461 

0.4950 

1.869 

3950 

34.9 

0.2602 

0.4983 

1.663 

4000 

39.1 

0.2758 

0.5018 

1.478 

4050 

43.6 

0.2931 

0.5057 

1.312 

4100 

48.5 

0.3125 

0.5098 

1.163 

4150 

53.8 

0.3348 

0.5144 

1.028 

4200 

59.6 

0.3606 

0.5195 

0.905 
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TABLE    2     (CONTINUED) 
CATHODE    SPOT     PARAMETERS     FOR     A     100     AMPERE     ARC    ON     A 

COPPER    CATHODE 

Cathode 

Electric 

Electron 

Ion 

Total 

Temperature 

Field 

Current 

Current 

Current 

(Kelvin) 

(10s    vVm) 

Density 

Density 

Density 

(10s    A/m2) 

(10s    A/m=) 

(10s    A/m  =  ) 

3500 

3.055 

0.625 

0.753 

1.378 

3550 

3.318 

0.841 

0.915 

1.756 

3600 

3.598 

1.124 

1.108 

2.232 

3650 

3.896 

1.495 

1.335 

2.830 

3700 

4.214 

1.980 

1.602 

3.582 

3750 

4.555 

2.610 

1.917 

4.527 

3800 

4.921 

3.430 

2.288 

5.717 

3850 

5.315 

4.493 

2.723 

7.216 

3900 

5.741 

5.874 

3.236 

9.110 

3950 

6.205 

7.670 

3.841 

11.51 

4000 

6.712 

10.01 

4.559 

14.57 

4050 

7.270 

13.07 

5.412 

18.49 

4100 

7.889 

17.11 

6.435 

23.54 

4150 

8.583 

22.47 

7.670 

30.14 

4200 

9.368 

29.69 

9.177 

38.87 

Ca  thode 

Ion/Total 

Coulomb 

Conduc- 

Erosion 

Temperature 

Current 

Logarithm 

tivity 

Rate 

(Kelvin) 

Density 
Ratio 

(mho/m) 

(10  =  ^ 
atoms/m^-s) 

3500 

.5462 

10.5 

3309 

3.226 

3550 

.5212 

10.5 

3489 

3.686 

3600 

.4963 

10.5 

3679 

4.193 

3650 

.4717 

10.4 

3879 

4.750 

3700 

.4474 

10.4 

4091 

5.360 

3750 

.4235 

10.4 

4315 

6.025 

3800 

.4001 

10.4 

4553 

6.748 

3850 

.3773 

10.3 

4808 

7.530 

3900 

.3552 

10.3 

5082 

8.372 

3950 

.3337 

10.3 

5378 

9.276 

4000 

.3129 

10.3 

5701 

10.24 

4050 

.2928 

10.3 

6057 

11.27 

4100 

.2733 

10.3 

6452 

12.36 

4150 

.2545 

10.3 

6898 

13.50 

4200 

.2361 

10.3 

7409 

14.69 
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TABLE     3 

CATHODE    SP 

150     AMPERE     ARC    ON     A 

COPPER    CATHODE 

Cathode 

Plasma 

Cathode 

Transition 

Total    Arc 

Temperature 

Temperature 

Potential 

Potential 

Potential 

(Kelvin) 

(eV) 

(Volts) 

(Volts) 

(Volts) 

3500 

1.48 

22.66 

10.79 

33.45 

3550 

1.53 

21.33 

11.41 

32.74 

3600 

1.58 

20.15 

12.03 

32.18 

3650 

1.64 

19.09 

12.65 

31.74 

3700 

1.69 

18.16 

13.25 

31.41 

3750 

1.75 

17.33 

13.84 

31.17 

3800 

1.81 

16.60 

14.41 

31.01 

3850 

1.88 

15.96 

14.93 

30.89 

3900 

1.94 

15.41 

15.40 

30.81 

3950 

2.02 

14.94 

15.79 

30.73 

4000 

2.09 

14.54 

16.08 

30.63 

4050 

2.18 

14.22 

16.25 

30.47 

4100 

2.27 

13.97 

16.26 

30.23 

4150 

2.37 

13.80 

16.06 

29.85 

4200 

2.49 

13.70 

15.59 

29.30 

Cathode 

Partial 

Degree    o-f 

Return 

Spot 

Temperature 

Pressure 

Ionization 

Coe-f -fie  lent 

R  a  d  iu  s 

(Kelvin) 

(atm.) 

(10-*    m) 

3500 

10.9 

0.1589 

0.4716 

5.948 

3550 

12.6 

0.1674 

0.4742 

5.269 

3600 

14.5 

0.1763 

0.4768 

4.674 

3650 

16.6 

0.1855 

0.4794 

4.151 

3700 

19.0 

0.1952 

0.4821 

3.690 

3750 

21.6 

0.2055 

0.4849 

3.282 

3800 

24.5 

0.2164 

0.4877 

2.921 

3850 

27.6 

0.2281 

0.4907 

2.600 

3900 

31.1 

0.2407 

0.4938 

2.315 

3950 

34.9 

0.2544 

0.4971 

2.060 

4000 

39.1 

0.2696 

0.5006 

1.831 

4050 

43.6 

0.2864 

0.5044 

1.626 

4100 

48.5 

0.3054 

0.5084 

1.442 

4150 

53.8 

0.3271 

0.5129 

1.275 

4200 

59.6 

0.3522 

0.5180 

1.123 
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TABLE    3     (CONTINUED) 
CATHODE    SPOT    PARAMETERS    FOR    A     150    AMPERE    ARC    ON    A 

COPPER    CATHODE 

Cathode 

Temperature 

(Kelvin) 

Electric 

Field 
(10s    vVm) 

Electron 
Current 
Density 

(10e    A/m=) 

Ion 

Current 

Density 

(10s    A/m  =  ) 

Total 

Current 

Density 

(10s    A/m2) 

3500 
3550 
3600 
3650 
3700 
3750 
3800 
3850 
3900 
3950 
4000 
4050 
4100 
4150 
4200 

3.011 
3.270 
3.546 
3.840 
4.153 
4.489 
4.850 
5.238 
5.658 
6.115 
6.614 
7.163 
7.772 
8.454 
9.225 

0.615 
0.827 
1.105 
1.469 
1.944 
2.562 
3.364 
4.405 
5.755 
7.509 
9.792 
12.78 
16.71 
21.92 
28.92 

0.735 
0.893 
1.081 
1.303 
1.564 
1.871 
2.232 
2.657 
3.157 
3.747 
4.445 
5.277 
6.272 
7.473 
8.937 

1.349 
1.720 
2.185 
2.771 
3.508 
4.433 
5.596 
7.061 
8.912 
11.26 
14.24 
18.06 
22.98 
29.39 
37.86 

Cathode 

Temperature 

(Kelvin) 

Ion/Total 

Current 

Density 

Ratio 

Coulomb 
Logarithm 

Conduc- 
tivity 
(mho/m) 

Erosion 

Rate 

(lO27 

atoms/mx-s) 

3500 
3550 
3600 
3650 
3700 
3750 
3800 
3850 
3900 
3950 
4000 
4050 
4100 
4150 
4200 

0.5443 
0.5194 
0.4946 
0.4700 
0.4458 
0.4221 
0.3989 
0.3762 
0.3542 
0.3329 
0.3122 
0.2923 
0.2729 
0.2542 
0.2361 

10.5 
10.5 
10.5 
10.4 
10.4 
10.4 
10.4 
10.3 
10.3 
10.3 
10.3 
10.3 
10.3 
10.3 
10.3 

3279 
3458 
3645 
3842 
4051 
4272 
4506 
4757 
5027 
5318 
5636 
5985 
6372 
6808 
7307 

3.232 
3.692 
4.200 
4.759 
5.371 
6.038 
6.762 
7.547 
8.392 
9.299 
10.27 
11.30 
12.39 
13.54 
14.74 
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TABLE     4 

CATHODE    SP( 

175     AMPERE    ARC    ON    A 

COPPER    CATHODE 

Cathode 

Plasma 

Cathode 

Transition 

Total    Arc 

Temperature 

Temperature 

Potential 

Potential 

Potential 

(Kelvin) 

<eV) 

(Volts) 

(Volts) 

(Volts) 

3500 

1.48 

22.61 

11.69 

34.29 

3550 

1.53 

21.26 

12.37 

33.63 

3600 

1.58 

20.08 

13.04 

33.13 

3650 

1.63 

19.04 

13.72 

32.76 

3700 

1.69 

18.10 

14.39 

32.49 

3750 

1.75 

17.28 

15.04 

32.32 

3800 

1.81 

16.55 

15.66 

32.21 

3850 

1.87 

15.92 

16.24 

32.16 

3900 

1.94 

15.37 

16.76 

32.13 

3950 

2.01 

14.90 

17.21 

32.11 

4000 

2.09 

14.50 

17.55 

32.05 

4050 

2.17 

14.18 

17.76 

31.95 

4100 

2.26 

13.93 

17.80 

31.74 

4150 

2.37 

13.76 

17.63 

31.39 

4200 

2.48 

13.66 

17.18 

30.84 

Cathode 

Partial 

Degree    o-f 

Return 

Spot 

Temperature 

Pressure 

Ionization 

Coe-f-ficient 

Radius 

(Kelvin) 

(atm.) 

(10-*    m) 

3500 

10.9 

0.1578 

0.4714 

6.452 

3550 

12.6 

0.1662 

0.4739 

5.711 

3600 

14.5 

0.1750 

0.4765 

5.066 

3650 

16.6 

0.1842 

0.4791 

4.499 

3700 

19.0 

0.1938 

0.4818 

3.999 

3750 

21.6 

0.2040 

0.4845 

3.556 

3800 

24.5 

0.2148 

0.4874 

3.166 

3850 

27.6 

0.2264 

0.4903 

2.819 

3900 

31.1 

0.2389 

0.4934 

2.509 

3950 

34.9 

0.2525 

0.4967 

2.233 

4000 

39.1 

0.2676 

0.5002 

1.986 

4050 

43.6 

0.2843 

0.5039 

1.763 

4100 

48.5 

0.3031 

0.5080 

1.563 

4150 

53.8 

0.3246 

0.5125 

1.382 

4200 

59.6 

0.3495 

0.5175 

1.218 
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TABLE    4     (CONTINUED) 
CATHODE    SPOT     PARAMETERS     FOR    A     175     AMPERE     ARC    ON    A 

COPPER    CATHODE 

Cathode 

Electric 

Electron 

Ion 

Total 

Temperature 

Field 

Current 

Current 

Current 

(Kelvin) 

(10e    vVm) 

Density 

Density 

Density 

(10e    A/m=) 

(10s    A/m  =  ) 

(10e    A/m  =  ) 

3500 

2.986 

0.609 

0.729 

1.338 

3550 

3.254 

0.822 

0.886 

1.708 

3600 

3.528 

1.098 

1.072 

2.171 

3650 

3.821 

1.460 

1.292 

2.752 

3700 

4.133 

1.932 

1.551 

3.483 

3750 

4.468 

2.546 

1.856 

4.402 

3800 

4.827 

3.342 

2.214 

5.557 

3850 

5.213 

4.376 

2.635 

7.011 

3900 

5.631 

5.716 

3.131 

8.847 

3950 

6.085 

7.456 

3.716 

11.17 

4000 

6.581 

9.722 

4.408 

14.13 

4050 

7.128 

12.68 

5.232 

17.91 

4100 

7.734 

16.58 

6.218 

22.79 

4150 

8.412 

21.74 

7.408 

29.15 

4200 

9.178 
Ion/Total 

28.67 

8.859 

37.53 

Cathode 

Coulomb 

Conduc- 

Erosion 

Temperature 

Current 

Logarithm 

tivity 

Rate 

(Kelvin) 

Density 
Ratio 

(mho/m) 

(10  =  TT 
a toms/m=-s) 

3500 

0.5448 

10.5 

3270 

3.233 

3550 

0.5188 

10.5 

3447 

3.694 

3600 

0.4940 

10.5 

3634 

4.203 

3650 

0.4695 

10.4 

3830 

4.762 

3700 

0.4453 

10.4 

4038 

5.374 

3750 

0.4216 

10.4 

4258 

6.042 

3800 

0.3984 

10.4 

4491 

6.767 

3850 

0.3758 

10.3 

4741 

7.552 

3900 

0.3539 

10.3 

5009 

8.398 

3950 

0.3326 

10.3 

5299 

9.306 

4000 

0.3120 

10.3 

5614 

10.28 

4050 

0.2921 

10.3 

5961 

11.31 

4100 

0.2728 

10.3 

6346 

12.40 

4150 

0.2542 

10.3 

6779 

13.55 

4200 

0.2361 

10.3 

7274 

14.75 
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TABLE    5 

CATHODE    SP< 

200    AMPERE    ARC    ON    A 

COPPER    CATHODE 

Cathode 

Plasma 

Cathode 

Transition 

Total    Arc 

Temperature 

Temperature 

Potential 

Potential 

Potential 

(Kelvin) 

<eV) 

(Volts) 

(Volts) 

(Volts) 

3500 

1.48 

22.48 

12.53 

35.01 

3550 

1.53 

21.21 

13.26 

34.47 

3600 

1.58 

20.03 

13.99 

34.02 

3650 

1.63 

18.99 

14.72 

33.71 

3700 

1.69 

18.06 

15.44 

33.50 

3750 

1.74 

17.24 

16.15 

33.38 

3800 

1.80 

16.51 

16.82 

33.34 

3850 

1.87 

15.88 

17.46 

33.34 

3900 

1.94 

15.33 

18.03 

33.37 

3950 

2.01 

14.86 

18.53 

33.39 

4000 

2.08 

14.47 

18.92 

33.39 

4050 

2.17 

14.15 

19.17 

33.32 

4100 

2.26 

13.90 

19.24 

33.15 

4150 

2.36 

13.73 

19.09 

32.82 

4200 
Cathode 

2.47 

13.63 

18.65 

32.28 

Partial 

Degree    o-f 

Return 

Spot 

Temperature 

Pressure 

Ionization 

Coe-f-ficient 

Radius 

(Kelvin) 

(atm.) 

(10~4    m) 

3500 

10.9 

0.1575 

0.4713 

6.899 

3550 

12.6 

0.1652 

0.4736 

6.122 

3600 

14.5 

0.1739 

0.4762 

5.431 

3650 

16.6 

0.1830 

0.4788 

4.823 

3700 

19.0 

0.1926 

0.4815 

4.287 

3750 

21.6 

0.2027 

0.4842 

3.814 

3800 

24.5 

0.2135 

0.4871 

3.395 

3850 

27.6 

0.2250 

0.4900 

3.022 

3900 

31.1 

0.2374 

0.4931 

2.690 

3950 

34.9 

0.2510 

0.4964 

2.394 

4000 

39.1 

0.2659 

0.4998 

2.129 

4050 

43.6 

0.2825 

0.5036 

1.891 

4100 

48.5 

0.3012 

0.5076 

1.677 

4150 

53.8 

0.3226 

0.5121 

1.483 

4200 

59.6 

0.3473 

0.5171 

1.307 
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TABLE    5     (CONTINUED) 
CATHODE    SPOT    PARAMETERS     FOR    A    200    AMPERE    ARC    ON    A 

COPPER    CATHODE 

Cathode 

Electric 

Electron 

Ion 

Total 

Temperature 

Field 

Current 

Current 

Current 

(Kelvin) 

(10e    V/m) 

Density 
(10e    A/m  =  ) 

Density 
(10e    A/m  =  ) 

Density 
(10e    A/m  =  ) 

3500 

2.990 

0.610 

0.727 

1.338 

3550 

3.241 

0.818 

0.880 

1.698 

3600 

3.515 

1.093 

1.065 

2.159 

3650 

3.806 

1.453 

1.284 

2.737 

3700 

4.117 

1.923 

1.541 

3.464 

3750 

4.450 

2.533 

1.844 

4.377 

3800 

4.808 

3.326 

2.199 

5.525 

3850 

5.193 

4.353 

2.617 

6.970 

3900 

5.609 

5.685 

3.110 

8.795 

3950 

6.061 

7.414 

3.691 

11.11 

4000 

6.556 

9.665 

4.379 

14.04 

4050 

7.100 

12.61 

5.197 

17.80 

4100 

7.703 

16.47 

6.176 

22.65 

4150 

8.378 

21.60 

7.357 

28.95 

4200 

9.141 

28.47 

8.797 

37.27 

Cathode 

Ion/Total 

Coulomb 

Conduc- 

Erosion 

Temperature 

Current 

Logarithm 

tivity 

Rate 

(Kelvin) 

Density 
Ratio 

(mho/m) 

(10=^ 

atoms/m=-s) 

3500 

0.5438 

10.5 

3267 

3.234 

3550 

0.5183 

10.5 

3439 

3.696 

3600 

0.4935 

10.5 

3625 

4.205 

3650 

0.4690 

10.4 

3820 

4.765 

3700 

0.4449 

10.4 

4027 

5.377 

3750 

0.4212 

10.4 

4246 

6.045 

3800 

0.3981 

10.4 

4479 

6.771 

3850 

0.3755 

10.3 

4727 

7.556 

3900 

0.3536 

10.3 

4994 

8.403 

3950 

0.3323 

10.3 

5283 

9.312 

4000 

0.3118 

10.3 

5597 

10.28 

4050 

0.2919 

10.3 

5942 

11.32 

4100 

0.2727 

10.3 

6325 

12.41 

4150 

0.2541 

10.3 

6755 

13.56 

4200 

0.2360 

10.3 

7247 

14.77 
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Figure    5.5      Potentials    for  a   200    Ampere    Arc 
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TABLE    6 

1 

CATHODE    SPi 

300    AMPERE    ARC    ON    A 

COPPER    CATHODE 

Cathode 

Plasma 

Cathode 

Transition 

Total    Arc 

Temperature 

Temperature 

Potential 

Potential 

Potential 

(Kelvin) 

(eV) 

(Volts) 

(Volts) 

(Volts) 

3500 

1.47 

22.32 

15.45 

37.77 

3550 

1.52 

21.06 

16.36 

37.42 

3600 

1.57 

19.90 

17.28 

37.17 

3650 

1.62 

18.86 

18.20 

37.06 

3700 

1.68 

17.94 

19.11 

37.05 

3750 

1.74 

17.13 

20.00 

37.13 

3800 

1.80 

16.41 

20.87 

37.28 

3850 

1.86 

15.78 

21.69 

37.47 

3900 

1.93 

15.24 

22.45 

37.69 

3950 

2.00 

14.77 

23.11 

37.89 

4000 

2.07 

14.38 

23.66 

38.04 

4050 

2.15 

14.07 

24.04 

38.11 

4100 

2.25 

13.82 

24.23 

38.05 

4150 

2.34 

13.65 

24.15 

37.80 

4200 

2.46 

13.55 

23.75 

37.30 

Cathode 

Partial 

Degree    o-f 

Return 

Spot 

Temperature 

Pressure 

Ionization 

Coe-f-ficient 

Radius 

(Kelvin) 

(atm.) 

(10-"*    m) 

3500 

10.9 

0.1550 

0.4706 

8.513 

3550 

12.6 

0.1626 

0.4730 

7.554 

3600 

14.5 

0.1712 

0.4755 

6.701 

3650 

16.6 

0.1801 

0.4781 

5.951 

3700 

19.0 

0.1896 

0.4808 

5.290 

3750 

21.6 

0.1995 

0.4835 

4.707 

3800 

24.5 

0.2101 

0.4863 

4.189 

3850 

27.6 

0.2214 

0.4892 

3.730 

3900 

31.1 

0.2336 

0.4923 

3.321 

3950 

34.9 

0.2470 

0.4955 

2.956 

4000 

39.1 

0.2616 

0.4990 

2.629 

4050 

43.6 

0.2779 

0.5027 

2.336 

4100 

48.5 

0.2963 

0.5067 

2.071 

4150 

53.8 

0.3172 

0.5111 

1.832 

4200 

59.6 

0.3415 

0.5160 

1.616 
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TABLE    6     (CONTINUED) 
CATHODE    SPOT    PARAMETERS    FOR    A    300    AMPERE    ARC    ON    A 

COPPER    CATHODE 

Cathode 

Electric 

Electron 

Ion 

Total 

Temperature 

Field 

Current 

Current 

Current 

(Kelvin) 

(10s    V/m) 

Density 

Density 

Density 

(10s    A/m=> 

(10s    A/m  =  ) 

(10s    A/m=) 

3500 

2.959 

0.603 

0.715 

1.318 

3550 

3.207 

0.808 

0.865 

1.673 

3600 

3.478 

1.080 

1.047 

2.127 

3650 

3.767 

1.435 

1.261 

2.696 

3700 

4.075 

1.898 

1.514 

3.412 

3750 

4.404 

2.500 

1.811 

4.311 

3800 

4.758 

3.280 

2.161 

5.441 

3850 

5.139 

4.292 

2.572 

6.863 

3900 

5.551 

5.603 

3.055 

8.658 

3950 

5.998 

7.304 

3.625 

10.93 

4000 

6.487 

9.516 

4.300 

13.82 

4050 

7.025 

12.40 

5.103 

17.51 

4100 

7.621 

16.20 

6.063 

22.26 

4150 

8.287 

21.22 

7.221 

28.44 

4200 

9.040 

27.95 

8.632 

36.58 

Cathode 

Ion/Total 

Coulomb 

Conduc- 

Erosion 

Temperature 

Current 

Logarithm 

tivity 

Rate 

(Kelvin) 

Density 
Ratio 

(mho/m) 

QO27 
atoms/m^-s) 

3500 

0.5423 

10.5 

3246 

3.238 

3550 

0.5169 

10.5 

3416 

3.701 

3600 

0.4922 

10.5 

3600 

4.211 

3650 

0.4678 

10.4 

3794 

4.771 

3700 

0.4438 

10.4 

3999 

5.385 

3750 

0.4202 

10.4 

4216 

6.054 

3800 

0.3971 

10.4 

4446 

6.781 

3850 

0.3747 

10.3 

4692 

7.568 

3900 

0.3528 

10.3 

4956 

8.417 

3950 

0.3317 

10.3 

5241 

9.328 

4000 

0.3113 

10.3 

5551 

10.30 

4050 

0.2915 

10.3 

5891 

11.34 

4100 

0.2724 

10.3 

6269 

12.44 

4150 

0.2539 

10.3 

6693 

13.59 

4200 

0.2360 

10.3 

7176 

14.80 
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Figure    5.6      Potentials    for   a    300    Ampere    Arc 
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TABLE     7 

CATHODE    SPi 

400     AMPERE     ARC    ON    A 

COPPER    CATHODE 

Cathode 

Plasma 

Cathode 

Transition 

Total    Arc 

Temperature 

Temperature 

Potential 

Potential 

Potential 

(Kelvin) 

(eV) 

(Volts) 

(Volts) 

(Volts) 

3500 

1.47 

22.22 

17.91 

40.13 

3550 

1.52 

20.97 

18.97 

39.94 

3600 

1.57 

19.81 

20.05 

39.86 

3650 

1.62 

18.78 

21.12 

39.91 

3700 

1.67 

17.87 

22.20 

40.07 

3750 

1.73 

17.06 

23.25 

40.31 

3800 

1.79 

16.35 

24.28 

40.63 

3850 

1.85 

15.72 

25.25 

40.98 

3900 

1.92 

15.18 

26.16 

41.34 

3950 

1.99 

14.72 

26.97 

41.69 

4000 

2.07 

14.33 

27.64 

41.98 

4050 

2.15 

14.02 

28.15 

42.16 

4100 

2.24 

13.77 

28.42 

42.20 

4150 

2.34 

13.60 

28.41 

42.01 

4200 

2.45 

13.50 

28.03 

41.53 

Cathode 

Partial 

Degree    o-f 

Return 

Spot 

Temperature 

Pressure 

Ionization 

Coe-f-ficient 

R  a  d  iu  s 

(Kelvin) 

(atm.) 

(10-4    m) 

3500 

10.9 

0.1535 

0.4702 

9.874 

3550 

12.6 

0.1610 

0.4726 

8.762 

3600 

14.5 

0.1695 

0.4751 

7.772 

3650 

16.6 

0.1784 

0.4777 

6.903 

3700 

19.0 

0.1877 

0.4803 

6.136 

3750 

21.6 

0.1976 

0.4830 

5.459 

3800 

24.5 

0.2081 

0.4858 

4.860 

3850 

27.6 

0.2193 

0.4888 

4.327 

3900 

31.1 

0.2314 

0.4918 

3.853 

3950 

34.9 

0.2446 

0.4950 

3.430 

4000 

39.1 

0.2591 

0.4984 

3.051 

4050 

43.6 

0.2752 

0.5021 

2.710 

4100 

48.5 

0.2934 

0.5061 

2.404 

4150 

53.8 

0.3141 

0.5105 

2.127 

4200 

59.6 

0.3380 

0.5154 

1.876 
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TABLE    7     (CONTINUED) 
CATHODE    SPOT    PARAMETERS    FOR    A    400    AMPERE    ARC    ON    A 

COPPER    CATHODE 

Cathode 

Electric 

Electron 

Ion 

Total 

Temperature 

Field 

Current 

Current 

Current 

(Kelvin) 

(10s    V/m) 

Density 
(10e    A/m=) 

Density 
(10e    A/m2) 

Density 
(10e    A/m2) 

3500 

2.940 

0.599 

0.707 

1.306 

3550 

3.187 

0.803 

0.856 

1.659 

3600 

3.456 

1.072 

1.036 

2.108 

3650 

3.743 

1.424 

1.248 

2.672 

3700 

4.049 

1.883 

1.498 

3.3822 

3750 

4.377 

2.480 

1.792 

4.272 

3800 

4.728 

3.253 

2.138 

5.391 

3850 

5.107 

4.256 

2.544 

6.800 

3900 

5.516 

5.554 

3.022 

8.577 

3950 

5.961 

7.238 

3.586 

10.83 

4000 

6.446 

9.427 

4.254 

13.68 

4050 

6.980 

12.28 

5.047 

17.33 

4100 

7.572 

16.03 

5.996 

22.03 

4150 

8.234 

21.00 

7.141 

28.14 

4200 

8.981 

27.64 

8.534 

36.17 

Cathode 

Ion/Total 

Coulomb 

Conduc- 

Erosion 

Temperature 

Current 

Logarithm 

tivity 

Rate 

(Kelvin) 

Density 
Ratio 

(mho/m) 

(10=^ 

atoms/m=-s) 

3500 

0.5415 

10.5 

3234 

3.240 

3550 

0.5161 

10.5 

3403 

3.704 

3600 

0.4915 

10.5 

3586 

4.214 

3650 

0.4671 

10.4 

3779 

4.775 

3700 

0.4431 

10.4 

3983 

5.389 

3750 

0.4196 

10.4 

4198 

6.059 

3800 

0.3966 

10.4 

4427 

6.787 

3850 

0.3742 

10.3 

4671 

7.575 

3900 

0.3524 

10.3 

4933 

8.425 

3950 

0.3313 

10.3 

5216 

9.337 

4000 

0.3109 

10.3 

5524 

10.31 

4050 

0.2912 

10.3 

5861 

11.35 

4100 

0.2722 

10.3 

6236 

12.45 

4150 

0.2538 

10.3 

6656 

13.61 

4200 

0.2359 

10.3 

7135 

14.82 
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TABLE    8 
CATHODE    SPOT    PARAMETERS    FOR    STABLE    ARCS 

Arc   Current 

Parameter 

200 

300 

400 

Cathode    Temperature 
(Kelvin) 

3850 

3700 

3600 

Plasma   Temperature   (eV) 

L87 

L68 

1.57 

Cathode    Potential  (Volts) 

15.88 

17.94 

19.81 

Spot   Radius    (10-4    m) 

3.022 

5.290 

7.772 

Ion   Current   Density   (10 8 
A/m2) 

2.617 

1.514 

1.036 

Electron  Current   Density 
(108    A/m2) 

4.353 

L898 

1.072 

Total  Current   Density   (108 
A/m2) 

6.970 

3.412 

2.108 

Erosion  Rate   (102  7 
atoms /m2 -s) 

7.556 

5.385 

4.214 
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VI.        CQHCLPSIQHS    AfiJD.   RECOMMENDATIONS 

This  study  has  attempted  to  answer  some  of  the  questions 
associated  with  the  cathode  phenomena  in  a  unipolar  arc. 
However,  a  more  detailed  analysis  of  the  cathode  spot  must  be 
accomplished  before  a  consistent  understanding  of  the  diverse 
phenomena  occurring  within  the  arc  can  be  achieved.  The 
stationary  model  presented  her  is  only  a  beginning  and  a  tool 
for  the  development  of  future  models.  It  can  not  answer  the 
questions  regarding  arc  initiation,  extinction,  and  spot  motion 
that  will  have  to  be  addressed  in  future  models  if  they  are  to 
be  accepted.  Additional  research  with  stationary  models  should 
consider  including  the  reverse  electron  current  from  the  plasma 
back  to  the  metal,  changing  from  a  constant  cathode  temperature 
in  the  arc  radius  to  a  constant  heat  flux  with  a  varying 
temperature,  considering  the  ohmic  heating  of  the  cathode  spot 
by  the  electron  current,  and  varying  the  cathode  material 
parameters  to  determine  their  influence  on  the  arc  stability.  It 
is  hoped  that  this  study  will  help  to  stimulate  further  interest 
and    research   into    the    cathode    spot  of  a    unipolar   arc. 
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APPENDIX    A       GEHERAf,    EMISSION    PROGRAM 

****************************************** 
* 

EMM1.EXE 


PURPOSE 

CALCULATE    THE    ELECTRON   CURRENT    DENSITY 
EMITTED   FOR    A   RANGE    OF    ELECTRIC    FIELDS    AT 
THE    SURFACE    OF    A   CATHODE    WITH   A   GIVEN    WORK 
FUNCTION    AND   CATHODE   TEMPERATURE. 

DESCRIPTION   OF    VARIABLES 


T 

F 

FO 

PHI 

DF 

J 

ECODE 


TYP 


* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

REAL*8    T.F.FO.PHLDF.J 

INTEGER    I,N,ECODE,TYP 

OPEN(4,FILE='CON') 

OPEN(6,FrLE='EMMLDAT',STATUS  =  'NEW') 

WRITEC*,*)   'ENTER    THE    NUMBER    OF    DATA   POINTS:    ' 

READ(*,*)   N 

WRITE(*,*)  'ENTER  THE  ABSOLUTE  TEMPERATURE  (K):  ' 

READ(*,*)  T 

WRITE(*,*)   'ENTER    INITIAL    ELECTRIC    FIELD   (V/M):    ' 

READ(*,*)   FO 

WRITE(*,*)   'ENTER    ELECTRIC    FIELD   INCREMENT    (V/M):    ' 

READ(*,*)  DF 


CATHODE  TEMPERATURE 

ELECTRIC  FIELD 

INITIAL  ELECTRIC  FIELD 

WORK  FUNCTION 

ELECTRIC  FIELD  INCREMENT 

CURRENT  DENSITY 

ERROR  CODE 

0  =  GOOD 

1  =  ERROR 
TYPE  OF  EMISSION 

1=  THERMIONIC 

2-  INTERMEDIATE 

3-  FIELD 
4=    NONE 


SUBROUTINES    REQUIRED 
EMMAPP(T,F,PHI,TYP) 
TFEMMl(T,F,PHLJ,ECODE) 
IFTEMMl(T,F,PHI,J,ECODE) 
FTEMMl(T,F,PHLJ,ECODE) 
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WRITE(*,*)  'ENTER  WORK  FUNCTION  (eV):  ' 

READ(*i*)  PHI 

F=F0 

J=0. 

DO  20  1=1, N 

CALL  EMMAPP(T,F,PHI,TYP) 
IF(TYP.EQ.l)  CALL  TFEMM1(T,F,PHI, J.ECODE) 
IFCTYP.EQ.2)  CALL  IFTEMM1(T,F,PHI, J.ECODE) 
IF(TYP.EQ.3)  CALL  FTEMM1(T,F,PHI, J.ECODE) 
IFCTYP.EQ.4)  GOTO  10 
WRITE(6,900)  F,  J,  ECODE,  TYP 
10     F=F  +  DF 
20  CONTINUE 

STOP 
900  FORMAT(lX,D12.6,5X,D12.6(5X,I5,5X,I5) 

END 


APPENDIX    B       GENERAL    EMISSION    SOBROPTINR 

\^  ^  *^  *^  -^  ^S  ■^  ^  ^^  ^v  J(^  >^  ^^  ^^  ^  ^  -^  ^  ^  ^  'T'   ^  ^  ^  ^1^  ^r-  'T'  nr-   ■'T-  ^  ^S  HP   T^  *P  O^  t^  *t^  m^  ^  ^  *  ^P  t^  ^  *P  O^  T^  'K  5|>  3|C 

C  * 

C  *  SUBROUTINE    EMM 

C  * 

C  *  PURPOSE 

C  *  CALCULATE    THE    ELECTRON   CURRENT    DENSITY 

C  *  FOR    A   GIVEN   ELECTRIC    FIELD   AT   THE   SURFACE 

C  *  OF    A   CATHODE    WITH   A   GIVEN    WORK   FUNCTION    AND 

C  *  CATHODE    TEMPERATURE. 

C  * 

C  *  DESCRIPTION   OF    VARIABLES 

C  *  T                     CATHODE    TEMPERATURE 

C  *  F                     ELECTRIC    FIELD 

C  *  PHI                WORK   FUNCTION 

C  *  J                     CURRENT    DENSITY 

C  *  ECODE           ERROR    CODE 

C  *  0    =    GOOD 

C  *  1  =    ERROR 

C  *  TYP                TYPE    OF    EMISSION 

C  *  1=    THERMIONIC 

C  *  2-    INTERMEDIATE 

C  *  3=    FIELD 

C  *  4=    NONE 

C  * 

C  *  SUBROUTINES    REQUIRED 

C  *  EMMAPP(T,F,PHI,TYP) 

C  *  TFEMMl(T,F,PHLJ,ECODE) 

C  *  IFTEMM1(T,F,PHLJ,EC0DE) 

C  *  FTEMMl(T,F,PHI,J,ECODE) 

C  * 

C 

SUBROUTINE    EMM(T,F,PHI, J, ECODE, TYP) 

REAL*8    T.F.PHLJ 

INTEGER    ECODE, TYP 

CALL    EMMAPP(T,F,PHI,TYP) 

IF(TYP.EQ.l)  CALL   TFEMM1(T,F,PHI, J,ECODE) 

IF(TYP.EQ.2)  CALL    IFTEMM1(T,F,PHI, J.ECODE) 

IFCTYP.EQ.3)   CALL    FTEMMKT.F.PHI, J.ECODE) 

RETURN 

END 
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APPENDIX    C       KMTSSIQH    APPf.  THARTT.  TTY    RIIBROnTINK 


c 

************************************** 

c 

* 

c 

* 

SUBROUTINE 

EMMAPP 

c 

* 

c 

* 

PURPOSE 

c 

* 

DETERMINES    APPLICABILITY   OF   THE   SUBROUTINES 

c 

* 

TFEMM1,    : 

EFTEMM1,    AND   FTEMM1 

c 

* 

c 

* 

DESCRIPTION 

OF    VARIABLES 

c 

* 

T 

CATHODE   TEMPERATURE 

c 

* 

FO 

ELECTRIC   FIELD   AT   CATHODE   SURFACE 

c 

* 

PHIO 

WORK   FUNCTION 

c 

* 

TYP 

TYPE   OF    EMISSION 

c 

* 

1=    THERMIONIC 

c 

* 

2-    INTERMEDIATE 

c 

* 

3=    FIELD 

c 

* 

4=    NONE 

c 

* 

c 

* 

SUBROUTINES 

REQUIRED 

c 

* 

NONE 

c 

* 

c 
c 

X  .  j   .  i .-  .  i   j.-  .1   i,  x,  ^l.  ,x,  -x.-  x.  x.  .X/  %x-  .V  -X-  X-  •X''  -X*  •*•  -X  ~X-"  -X"  -X/  -X*'  "X-  X-  -X-  -X-  -X-  ■X'  'X-  'X-  -X"  'X-  ~X  X-  X/  "X-  -X/  X-  -X-  -X-  -X'  -X-  'X/  "X/  -X-  X/ 
^p.  ^i  ^i  ^^  ^^  ^i  ^^  ^^  ^^.  -^»  ^^  ^n  .^v  ^i  >^  ^p-  ^i  ^^  *^  *^  i^  *^>  <^  ^^  "X1  ^^  ^n  *^  ^n  ^^  *^*  ^H1  ^^  ^^  ^*  ^^  ^^  ^^  ^n  ^*  t*  ^n  ^^  ^^  ^^  ^n  *^  ^n  *T*  t* 

SUBROUTINE    EMMAPP(T,FO,PHIO,TYP) 
REAL*8    T,F,PHI,M,E,K,H,EPS,PLPHI0,D,D1 
REAL*8    F0,E1,F1,KT0,KT,Y,YE,TY,TE,V,VE,ETA,ETA0 
REAL*8    TEST1,TEST2,TEST3,TEST4 
INTEGER    ECODEl,ECODE2,ECODE3,TYP,N 
M=9.1095E-31 
E=1.6022E-19 
K=1.3807E-23 
H=1.0546E-34 
EPS  =  8.8542E-12 
PI=3. 141592654 
E1=(4*PI*EPS*H)**2/(M*E**4) 
F1=E1**2*E**3/(4*PI*EPS) 
F=F0*F1 
KT0=K*T 
KT=KT0*E1 
PHI=PHI0*E1*E 
N=0 

D=F**(3./4.)/(PI*KT) 
IF(l-D)   20,20,10 
10    TEST1=  DLOG((l- D)/D)-  1/(D*(1- D)) 

TEST2  =  -PI/(F**(3./4.))*(PHI-DSQRT(F)) 
TEST3=DLOG((l-D)/D)-l/(l-D) 
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TEST4=-PI/(F**(l./8.)) 

IF((TEST-1.GT.TEST2).AND.(TEST3.GT.TEST4))   THEN 

ECODE1=0 
ELSE 

EC0DE1=1 
ENDIF 
GOTO   30 
20    ECODEl-1 

30    ETA0  =  -F**2/(8*KT**2) 
Y=DSQRT(F)/PHI 
CALL    VTY(Y,V,TY,IER1) 
40    YE=DSQRT(F)/(-ETA0) 

CALL    VTY(YE,VE,TE,IER2) 
ETA=-F**2/(8*KT**2*TE**2) 
IF(DABS(ETA-ETAO).GT.    1.0D-6)   THEN 

ETA0  =  ETA-(ETA-ETA0)/2.1 

N=N+1 

IF(N.GT.200)   THEN 
ECODE2=l 
GOTO    50 

ENDIF 

GOTO    40 
ENDIF 

D=2*DSQRT(2/YE)*TE/PI 
TEST1=1/YE 

TEST2=1   +    F**(1./4.)*D/(PI*(D-1)) 
TEST3  =  -F**2/(8*KT**2*TE**2) 
TEST4  =  -PHI+KT/(1-F/(2*DSQRT(2*PHI)*KT*TY)) 
IF((TEST1.GT.TEST2).AND.(TEST3.GT.TEST4))   THEN 

ECODE2=0 
ELSE 

ECODE2-1 
ENDIF 
50    D=(2*DSQRT(2.0*PHI)*PI*KT)/F 
D1=4*DSQRT(2*PHI)*PHI/(3*F) 
TEST1=PHI-DSQRT(F) 

TEST2  =  F**(3./4.)/PI   +    KT/(1-D*TY/PI) 
TEST3=1-D*TY/PI 

TEST4  =  KT*DSQRT((V*PHI/F)*DSQRT(2*PHI)/(PHI**2-F)) 
IF((TESTLGT.TEST2).AND.(TEST3.GT.TEST4))   THEN 

ECODE3r0 
ELSE 

ECODE3=l 
ENDIF 
IF(ECODELEQ.O)   THEN 

TYP=1 
ELSE 

IFCECODE2.EQ.0)   THEN 
TYP-2 

ELSE 

IFCECODE3.EQ.0)   THEN 
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TYP=3 
ELSE 

TYP=4 
ENDIF 
ENDIF 
ENDIF 
RETURN 
END 
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APPENDIX    D       1^S_   EMISSION    SPBROOTINK 


c 

********************************************** 

c 

* 

c 

* 

SUBROUTINE 

TFEMM1 

c 

* 

c 

* 

PURPOSE 

c 

* 

CALCULATES  THE  ELECTRON  CURRENT  DENSITY 

c 

* 

FOR  THERMIONIC  FIELD  EMISSION 

c 

* 

c 

* 

DESCRIPTION 

OF  VARIABLES 

c 

* 

T 

CATHODE  TEMPERATURE 

c 

* 

FO 

ELECTRIC  FIELD  AT  CATHODE  SURFACE 

c 

* 

PHIO 

WORK  FUNCTION 

c 

* 

JO 

CURRENT  DENSITY 

c 

* 

ECODE 

ERROR  CODE 

c 

* 

0=  GOOD 

c 

* 

1=  ERROR 

c 

* 

c 

* 

SUBROUTINES 

REQUIRED 

c 

* 

NONE 

c 

* 

c 
c 

************************************************** 

SUBROUTINE    TFEMM1(T,F0,  PHIO,  JO,  ECODE) 
REAL*8    T,F,PHI,J,M,E,K,H,EPS,PLPHIO,D 
REAL*8    F0,J0,E1,F1,J1,KT0,KT 
REAL*8    TEST1,TEST2,TEST3,TEST4 
INTEGER    ECODE 
M=9.1095E-31 
E=L6022E-19 
K=L3807E-23 
H-1.0546E-34 
EPS  =  8.8542E-12 
PI=3. 141592654 
E1=(4*PI*EPS*H)**2/(M*E**4) 
F1=E1**2*E**3/(4*PI*EPS) 
J1=E1**2*H**3/(M*E) 
F=F0*F1 
KT0=K*T 
KT=KT0*E1 
PHI=PHI0*E1*E 
D=F**(3./4.)/(PI*KT) 

jr(KT**2/(2*PI**2))*((PI*D)/DSIN(PI*D))* 
1      DEXP(-(PHI-DSQRT(F))/KT) 
J0=J/J1 

IF(l-D)   20,20,10 
10    TESTIS DLOG((l- D)/D)-  1/(D*(1- D)) 
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TEST2  =  -PI/(F**(3./4.))*(PHI-DSQRT(F)) 
TEST3  =  DL0G((1-  D)/D)- 1/(1-  D) 
TEST4=-PI/(F**(l./8.)) 

IF((TEST1.GT.TEST2).AND.(TEST3.GT.TEST4))  THEN 

ECODE^O 
ELSE 

ECODE-1 
ENDIF 
RETURN 
20  ECODE=l 
RETURN 
END 
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APPENDIX    E       FIELD-THERMIONIC    EMISSION    SOBROPTINE 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

* 
* 

K*************************^ 

SUBROUTINE 

FTEMM1 

* 
* 

PURPOSE 

CALCULATES  THE  ELECTRON  CURRENT  DENSITY 
FOR  FIELD-THERMIONIC  EMISSION 

* 
* 

* 
* 

* 

* 

* 

DESCRIPTION 
T 
FO 

PHIO 
JO 
ECODE 

OF  VARIABLES 
CATHODE  TEMPERATURE 
ELECTRIC  FIELD  AT  CATHODE  SURFACE 
WORK  FUNCTION 
CURRENT  DENSITY 
ERROR  CODE 

0=  GOOD 

1=  ERROR 

c 
c 
c 
c 
c 

* 

SUBROUTINES  REQUIRED 
VTY(Y,V,TY,IER) 

SUBROUTINE    FTEMM1(T,F0, PHIO, JO, ECODE) 

REAL*8    T,F,PHI,J,M,E,K,H,EPS,PI,PHI0,D,D1,TY 

REAL*8    TEST1,TEST2,TEST3,TEST4,Y,V,C1,C2,C3 

REAL*8    F0,J0,E1,F1,J1,KT0,KT 

INTEGER    ECODE 

M=9.1095E-31 

E=1.6022E-19 

K=1.3807E-23 

H-1.0546E-34 

EPS  =  8.8542E-12 

PI=3. 141592654 

E1=(4*PI*EPS*H)**2/(M*E**4) 

F1=E1**2*E**3/(4*PI*EPS) 

J1=E1**2*H**3/(M*E) 

F=F0*F1 

KT0=K*T 

KT=KT0*E1 

PHI=PHI0*E1*E 

D=(2*DSQRT(2.0*PHI)*PI*KT)/F 

D1=4*DSQRT(2*PHI)*PHI/(3*F) 

Y=DSQRT(F)/PHI 

CALL    VTY(Y,V,TY,IER) 

J=F**2*(D*TY/DSIN(D*TY))*DEXP(-D1*V)/(16*PI**2*PHI*TY**2) 

J0=J/J1 
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TEST1=PHI-DSQRT(F) 

TEST2  =  F**(3./4.)/PI  +  KT/(1-D*TY/PI) 

TEST3=1-D*TY/PI 

TEST4  =  KT*DSQRT((V*PHI/F)*DSQRT(2*PHI)/(PHI**2-F)) 

IF((TEST1.GT.TEST2).AND.(TEST3.GT.TEST4))  THEN 

ECODE=0 
ELSE 

ECODE^l 
ENDIF 
RETURN 
END 
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APPENDIX    F       COMPLETE    ELLIPTIC     INTEGRAL    Q£_  T_flK   l&X 

KI1LD. 

C 

*************************************** 

c 

C  SUBROUTINE  CEL1 

C 

C  PURPOSE 

C  CALCULATE  COMPLETE  ELLIPTIC  INTEGRAL  OF  FIRST 

C  KIND 

C 

C  USAGE 

C  CALL  CEL1(RES,AK,IER) 

C 

C  DESCRIPTION  OF  PARAMETERS 

C  RES    -  RESULT  VALUE 

C  AK     -  MODULUS  (INPUT) 

C  IER    -  RESULTANT  ERROR  CODE  WHERE 

C  IERrO   NO  ERROR 

C  IER=1   AK  NOT  IN  RANGE  -1  TO  +1 

C 

C  REMARKS 

C  THE  RESULT  IS  SET  TO  1.E75  IF  ABS(AK)  GE  1 

C  FOR  MODULUS  AK  AND  COMPLEMENTARY  MODULUS  CK, 

C  EQUATION  AK*AK+CK*CK=LO  IS  USED. 

C  AK  MUST  BE  IN  THE  RANGE  -1  TO  +1 

C 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

C  NONE 

C 

C  METHOD 

C  DEFINITION 

C  CEL1(AK)=INTEGRAL(1/SQRT((1+T*T)*(1+(CK*T)**2)), 

C  SUMMED   OVER   T    FROM   0    TO   INFINITY). 

C  EQUIVALENT    ARE    THE    DEFINITIONS 

C  CELl(AK)=INTEGRAL(l/(COS(T)SQRT(l+(CK*TAN(T))**2)) 

C  SUMMED  OVER  T  FROM  0  TO  PI/2), 

C  CEL1(AK)=INTEGRAL(1/SQRT(1-(AK*SIN(T))**2), 

C  SUMMED  OVER  T  FROM  0  TO  PI/2),  WHERE 

C  K=SQRT(1.-CK*CK). 

C  EVALUATION 

C  LANDENS  TRANSFORMATION  IS  USED  FOR  CALCULATION. 

C  REFERENCE 

C  R.BULIRSCH,    'NUMERICAL    CALCULATION    OF    ELLIPTIC 

C  INTEGRALS    AND    ELLIPTIC    FUNCTIONS',    HANDBOOK   SERIES 

C  SPECIAL    FUNCTIONS,    NUMERISCHE    MATHEMATIK   VOL.    7, 

C  1965,    PP.    78-90. 
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c 

[^  ^  ^r  ^h  ^  ^  ^-^F  t^  *1^  ^  3|5  ?|t  ^  >t^  ^  o^  ^k  /^C  3fC  y^C  )JC  3JC  >^t  ^  ^  ^  ^.  f   f    i(  )fC  ?(C  ?(t  JH.  JfC  )fC  ?fC  ?fC  >fC  ^C  )lC  )(C  JfC  )jC  ?(C  )^C  3(C  ^  )^C  >fC 

c 

SUBROUTINE   CELHRES,  AK.IER) 

REAL*8    ARI,GEO,AK,RES,AARI 

IER=0 

ARI=2. 

GEO=(0.5-AK)+0.5 

GEO=GEO+GEO*AK 

RES  =  0.5 

IF(GEO)l,2,4 

1  IER  =  1 

2  RES=LE38 
RETURN 

3  GEO=GEO*AARI 

4  GEO=DSQRT(GEO) 
GEO=GEO+GEO 
AARI^ARI 
ARI=ARI+GEO 
RES=RES+RES 
IF(GEO/AARI-  0. 9999  )3 ,5,5 

5  RES  =  RES/ARI*6.283185E0 
RETURN 

END 


APPENDIX  G   COMPLETE  ELLIPTIC  INTEGRAL  Q£_  2HD_  KIND 

C 

********************************************* 

c 

C  SUBROUTINE  CEL2 

C 

C  PURPOSE 

COMPUTES    THE    GENERALIZED   COMPLETE    ELLIPTIC 

C  INTEGRAL    OF    SECOND   KIND. 
C 

C  USAGE 

C  CALL    CEL2(RES,AK,A,B,IER) 
C 

C  DESCRIPTION  OF  PARAMETERS 

C  RES    -  RESULT  VALUE 

C  AK     -  MODULUS  (INPUT) 

C  A        CONSTANT  TERM  IN  NUMERATOR 

C  B      -  FACTOR  OF  QUADRATIC  TERM  IN  NUMERATOR 

C  IER    -  RESULTANT  ERROR  CODE  WHERE 

C  IER=0   NO  ERROR 

C  IER=1   AK  NOT  IN  RANGE  -1  TO  +1 
C 

C  REMARKS 

C  FOR  ABS(AK)  GE  1  THE  RESULT  IS  SET  TO  1.E75  IF  B 

C  IS  POSITIVE,  TO  -1.E75  IF  B  IS  NEGATIVE. 

C  SPECIAL  CASES  ARE 

C  K(K)  OBTAINED  WITH  A  =  1,  B  =  1 

C  E(K)  OBTAINED  WITH  A  =  1,  B  =  CK*CK  WHERE  CK  IS 

C  COMPLEMENTARY  MODULUS. 

C  B(K)  OBTAINED  WITH  A  =  1,  B  =  0 

C  D(K)  OBTAINED  WITH  A  =  0,  B  =  1 

C  WHERE  K,  E,  B,  D  DEFINE  SPECIAL  CASES  OF  THE 

C  GENERALIZED  COMPLETE  ELLIPTIC  INTEGRAL  OF  SECOND 

C  KIND  IN  THE  USUAL  NOTATION,  AND  THE  ARGUMENT  K  OF 

C  THESE  FUNCTIONS  MEANS  THE  MODULUS. 
C 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

C  NONE 
C 

C  METHOD 

C  DEFINITION 

C  RES  =  INTEGRAL((A+B*T*T)/(SQRT((1+T*T) 

C  *(1+(CK*T)**2))*(1+T*T)) 

C  SUMMED  OVER  T  FROM  0  TO  INFINITY). 

C  EVALUATION 

C  LANDENS  TRANSFORMATION  IS  USED  FOR  CALCULATION. 

C  REFERENCE 
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C  R.BULIRSCH,    'NUMERICAL   CALCULATION   OF    ELLIPTIC 

C  INTEGRALS    AND   ELLIPTIC   FUNCTIONS',    HANDBOOK  SERIES 

C  SPECIAL    FUNCTIONS,    NUMERISCHE    MATHEMATIK   VOL.    7, 

C  1965,    PP.    78-90. 

C 

C 

SUBROUTINE    CEL2(RES,AK, A,B,IER) 

REAL*8    RES,AK,A,B,ARI,GEO,A1,BO,AARI 

IER=0 

ARI-2. 

GEO=(0.5-AK)+0.5 

GEO=GEO+GEO*AK 

RES  =  A 

A1=A+B 

B0=B+B 

IF(GEO)l,2,6 

1  IER=1 

2  IF(B)3,8,4 

3  RES=-LE38 
RETURN 

4  RES=1.E38 
RETURN 

5  GEO=GEO*AARI 

6  GEO=DSQRT(GEO) 
GEO-GEO+GEO 
AARI-ARI 
ARI-ARI+GEO 
B0=B0+RES*GEO 
RES=A1 
B0=B0+B0 
AlrBO/ARI+Al 
IF(GEO/AARI-0. 9999)5, 7,7 

7  RES=A1/ARI 
RES=RES+0.5707963E0*RES 

8  RETURN 
END 
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APPENDIX  H   FQWLBR-KORDHTKH  ELLIPTIC  FUNCTIONS 

C  * 

C  *  SUBROUTINE    VTY 

C  * 

C  *  PURPOSE 

C  *  CALCULATES    THE    FOWLER-NORDHIEM    ELLIPTIC 

C  *  FUNCTIONS 

C  * 

C  *  DESCRIPTION   OF    VARIABLES 

C  *  Y  FUNCTION    ARGUMENT 

C  *  V  FOWLER-NORDHIEM   ELLIPTIC    V 

C  *  TY  FOWLER-NORDHIEM   ELLIPTIC    T 

C  *  IER  ERROR   CODE 

C  *  0=    GOOD 

C  *  1=    ERROR 

C  * 

C  *  SUBROUTINES    REQUIRED 

C  *  CEL1(RES,AK,IER) 

C  *  CEL2(RES,AK,A,B,IER) 

C  * 

C 

SUBROUTINE    VTY(Y, V,TY,IER) 
REAL*8    Y,V,TY,C1,C2,C3 
INTEGER    IER,IER1,IER2 
IFCY.GT.    1.0)   THEN 

C1=DSQRT((Y-1)/(2*Y)) 

CALL    CEL1(C2,C1,IER1) 

CALL    CEL2(C3,C1,1.0D00,L0D00-C1*C1,IER2) 

V=-DSQRT(Y/2)*(-2*C3    +    (Y+1)*C2) 

TY-1/DSQRT(1+Y)*((1+Y)*C3    -    Y*C2) 
ELSE 

C1=DSQRT((1-Y)/(1+Y)) 

CALL   CEL1(C2,C1,IER1) 

CALL   CEL2(C3,C1,1.0D00,1.0D00-C1*C1,IER2) 

V=DSQRT(1+Y)*(C3-Y*C2) 

TY=1/DSQRT(1+Y)*((1+Y)*C3    -Y*C2) 
ENDIF 
IF((IERl.EQ.l).OR.(IER2.EQ.l))   THEN 

IER=1 
ELSE 

IER=0 
ENDIF 
RETURN 
END 
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APPENDIX    I        INTERMEDIATE    F-T    EMISSION    SPBROOTINE 


c 
c 
c 
c 
c 

■k   ~k   A   A-  "i?  "A'   &  A-    A   A  A'  A"  A"   A'    A   A   'At'   A"   A    A1   "A  "A'  'A  A'  'A'   A    AAA'    A    A    A    A    A  A"    A  A"   A  A"  A-  ■*•  •*-    A  >A  <A^  A-   ■4-  A-  >A  A' 
-T^  ^n  ■^  ^n  ^n  ■^  ^n  •^  'T1  *^  ^p  •'T*  ^  ^n  'r1  *T-  -T^  ^^  ^*  •T*  ^*  •"T*  ^n  ^p  y^  ^n  ^-  -^  •'T*  ^  -^  'T"  -t*  *X"  -^  ^n  ^p  ^S  ^p  ^  ^S  "T1  ^»  ^S  ^^  -^  t*  -t*  ^  ^ 

A' 

* 
* 
* 

SUBROUTINE 

IFTEMM1 

PURPOSE 

c 

* 

CALCULATES    THE    ELECTRON   CURRENT    DENSITY 

c 
c 
c 

* 
* 
* 

FOR    INTERMEDIATE    FIELD-THERMIONIC    EMISSION 

DESCRIPTION 

OF    VARIABLES 

c 

* 

T 

CATHODE   TEMPERATURE 

c 

* 

FO 

ELECTRIC    FIELD   AT   CATHODE    SURFACE 

c 

* 

PHIO 

WORK   FUNCTION 

c 

* 

JO 

CURRENT    DENSITY 

c 

* 

ECODE 

ERROR   CODE 

c 

* 

0=    GOOD 

c 

* 

1=    ERROR 

c 

* 

c 

* 

SUBROUTINES 

REQUIRED 

c 
c 
c 
c 

* 

VTY(Y,V,TY,IER) 

T* 

SUBROUTINE    IFTEMM1(T,F0,  PHIO,  JO,  ECODE) 
REAL*8    T,F,PHI,J,M,E,K,H,EPS,PI,PHIO,D,TE,TY 
REAL*8    TEST1,TEST2,TEST3,TEST4,Y,V,C1,C2,C3,VE 
REAL*8    F0,J0,E1,F1,J1,KT0,KT,ETA,ETA0,YE,THETA 
INTEGER    ECODE 
M=9.1095E-31 
E=1.6022E-19 
K=1.3807E-23 
H=L0546E-34 
EPS  =  8.8542E-12 
PI=3. 141592654 
E1=(4*PI*EPS*H)**2/(M*E**4) 
F1=E1**2*E**3/(4*PI*EPS) 
J1=E1**2*H**3/(M*E) 
F=F0*F1 
KT0=K*T 
KT=KT0*E1 
PHI=PHI0*E1*E 
ETA0  =  -F**2/(8*KT**2) 
Y-DSQRT(F)/PHI 
CALL    VTY(Y,V,TY,IER1) 
20    YE=DSQRT(F)/(-ETAO) 

CALL    VTY(YE,VE,TE,IER2) 
ETA=-F**2/(8*KT**2*TE**2) 
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IF(DABS(ETA-ETAO).GT.    1.0D-6)   THEN 

ETAO=ETA-(ETA-ETA0)/2.1 

GOTO   20 
ENDIF 

THETA=3/(TE**2)    -    2*VE/(TE**3) 

J=F*DSQRT(KT*TE/(2*PI))/(2*PD*DEXP(-PHI/KT    +    F**2* 
1      THETA/(24*KT**3)) 
J0=J/J1 

D=2*DSQRT(2/YE)*TE/PI 
TESTlrl/YE 

TEST2=1   +    F**(1./4.)*D/(PI*(D-1)) 
TEST3  =  -F**2/(8*KT**2*TE**2) 
TEST4r-PHI+KT/(l-F/(2*DSQRT(2*PHI)*KT*TY)) 
IF((TEST1.GT.TEST2).AND.(TEST3.GT.TEST4))   THEN 

ECODE=0 
ELSE 

ECODE=l 
ENDIF 
RETURN 
END 
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APPENDIX  J   AJLC  PROGRAM 

c  * 

C  *      ARC. EXE 

C  * 

C  *      PURPOSE 

C  *  DETERMINES  THE  CHARACTERISTIC  VALUES  FOR 

C  *  THE  CATHODE  SPOT  OF  AN  ARC  BURNING  ON  A 

C  *  COPPER  PLATE 

C  * 

C  *      DESCRIPTION  OF  VARIABLES 

C  *  LO  HEAT  OF  VAPORIZATION 

C  *  PO  PARTIAL  PRESSURE  CONSTANT 

C  *  MI  MASS  OF  COPPER  ATOM 

C  *  EI  IONIZATION  POTENTIAL  ENERGY 

C  *  PHIO  WORK  FUNCTION  (EV) 

C  *  BETA  FRACTION  ENERGY  REMOVED  WITH 

C  *  COULOMB  INTERACTION 

C  *  XI  THERMAL  CONDUCTIVITY 

C  *  EER  ERROSION  ENERGY 

C  *  Z  AVERAGE  ION  CHARGE 

C  *  NO  AVOGADRO'S  NUMBER 

C  *  M  ELECTRON  MASS 

C  *  E  ELECTRON  CHARGE 

C  *  K  BOLTZMAN'S  CONSTANT 

C  *  H  PLANCK'S  CONSTANT 

C  *  EPS  PERMITIVITY  OF  FREE  SPACE 

C  *  ATP  DEGREE  OF  IONIZATION 

C  *  ETA  VAPORIZATION  RETURN  COEFFICIENT 

C  *  F  ELECTRIC  FIELD  AT  CATHODE  SURFACE 

C  *  G  EROSION  RATE 

C  *  PHI  SCHOTTKY  CORRECTED  WORK  FUNCTION 

C  *  PI  PI 

C  *  PTK  VAPOR  PARTIAL  PRESSURE 

C  *  R  RADIUS  OF  CATHODE  SPOT 

C  *  SIG  PLASMA  CONDUCTIVITY 

C  *  TK  CATHODE  TEMPERATURE  (K) 

C  *  U  POTENTIAL  DROP  IN  QUASINEUTRAL 

C  *  PLASMA 

C  *  VK  CATHODE  POTENTIAL  DROP 

C  *  VT  TOTAL  POTENTIAL  DROP  ACROSS  ARC 

C  *  J  TOTAL  CURRENT  DENSITY 

C  *  JE  ELECTRON  CURRENT  DENSITY 

C  *  JI  ION  CURRENT  DENSITY 

C  *  KTK  CATHODE  TEMPERATURE  (JOULES) 

C  *  KTP  PLASMA  TEMPERATURE  (JOULES) 

C  *  I  ARC  CURRENT 
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c 

* 

ND 

NUMBER   OF    DATA   SETS 

c 

* 

•      DTK 

CATHODE    TEMPERATURE    INCREMENT 

c 

* 

COULOG 

COULOMB   LOGARITHM 

c 

* 

ECODE 

ERROR    CODE 

c 

* 

0=    GOOD 

c 

* 

1=    ERROR 

c 

* 

TYP 

TYPE    OF    EMISSION 

c 

* 

1=    THERMIONIC 

c 

* 

2-    INTERMEDIATE 

c 

* 

3=    FIELD 

c 

* 

4=    NONE 

c 

* 

c 

* 

SUBROUTINES 

REQUIRED 

c 

* 

EMM(T,F,PHIO,  J,  ECODE, TYPE) 

c 

* 

c 
c 

-^  ^n1  ^'  "T^  "T*  'T'  ■T"  ^h"  ^-  -^  ^*  -^S  ^*  ^n1  -^  ^*  ^p  ^n  ^p>  ^F-  ^p  ^S  ^p-  ^P  ^p-  ^p>  ^p  ^^  ^p^  t^  ^n  ^n  ^F  ^p>  ^F  ^^  t^  m^  ^F  '!*■  'v*  ^n  ^n  ^p  ^p  't^  ^t 

REAL*8    LO,PO,MI,EI,PHIO,BETA,XI,EER,Z,NO,M,E,K,H,EPS 

REAL*8    ATP,ETA,F,G,PHI,PI,PTK,R,STG,TK,U,VK,VK1,VK2,VT 
REAL*8    Sl,S2,S3,S4,S5,J,JE,JLKTK,KTP,I,DTK,COULOG 
INTEGER    ECODE, TYP,N,ND 
OPEN(6,FILE='ARC.DAT',STATUS  =  'NEW') 
OPEN(7,FILE='PRN') 

OF    DATA   SETS:    ' 


WRITE(*,*)   'ENTER 

READ(*,*)   ND 

WRITE(*,*)   'ENTER 

READ(*,*)   KTP 

WRITE(*,*)   'ENTER 

READ(*,*)   TK 

WRITE(*,*)   'ENTER 

READ(*,*)   DTK 

WRITEC*,*)   'ENTER 

READ(*,*)   I 

WRITE(6,900)   ND,I 

WRITE(7,900)   ND,I 

MI=63.54 

L0=300300. 

EI=L2377D-18 

PHIO-4.4 

H=40. 

BETA=3.21 

Z=L8 

N0  =  6.0220D23 

M=9.1095D-31 

E  =  L6022D-19 

K=L3807D-23 

H=6.6262D-34 

PI=3. 141592654 

EPS-8.8542D-12 

EER=L0/N0 

S1=-L7520D4 


NUMBER 


INITIAL    PLASMA   TEMPERATURE    (EV):    ' 
INITIAL   CATHODE   TEMPERATURE    (K):    ' 
CATHODE    TEMPERATURE    INCREMENT    (K): 
ARC    CURRENT    (A):    ' 
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S2=-L21 
S3  =  0  .- 
S4=13.21 

S5  =  133. 322368421 
VK=15. 

MI=MI*L6726D-27 
KTP=E*KTP 
DO   20    N=1,ND 
KTK=K*TK 
JE=0. 

PTK=(10**(S1/TK   +S2*DLOG10(TK)+S3*TK+S4))*S5 
10    ATP=1/DSQRT(L    +   PTK/KTP*(2*PI*H**2/(M*KTP))**(3./2.)* 
1  DEXP(EI/KTP)) 

ETA=(8./3.    +  ATP-DSQRT(5.*PI*KTK/(6.*KTP)))/ 
1  (8./3.    +  ATP+DSQRT(5.*PI/6.)) 

JI=E*PTK*ETA*ATP/DSQRT(2*PI*MI*KTK) 
F  =  DSQRT(JI*DSQRT(8*VK*MI/E)/EPS*(1-JE/JI*DSQRT(M/MI))) 
PHI=PHI0*E-DSQRT(E**3*F/(4*PI*EPS)) 
CALL    EMM(TK,F,PHI0,JE,ECODE,TYP) 
R=  DSQRT(I/(PI*(  JE+ JI))) 

VK1=((JI/JE)*(EI+2*KTP)+BETA*KTP*(L+JI/JE)-2*KTK)/E 
G=PTK*(1-ETA)/DSQRT(2*PI*MI*KTK) 

VK2  =  4*XI*TK/(PI*R*JI)+((JE/JI)*(PHI+2*KTK)+(G/JI)*EER*E 
1  -EI-2*KTP+PHI)/E 

IF(DABS(VK2-VK1).GT.    1.0D-1)   THEN 
IF(VK2.GT.VK1)   THEN 
VK^VKl 
IF(DABS(VK2-VK1).GT.    100)   THEN 

KTP=KTP+100.*K 
ELSE 

KTP=KTP+DABS(VK2-VK1)*K 
ENDIF 
ELSE 

VK=VK2 

IF(DABS(VK2-VK1).GT.    100)   THEN 

KTP=KTP-10L*K 
ELSE 

KTP=  KTP-  DABS(VK2  -  VK1)*K 
ENDIF 
ENDIF 

IF(VK.LT.O)   VK^DABS(VK) 
GOTO   10 
ENDIF 

VK=(VKl+VK2)/2 
J=JE+JI 

COULOG=23.0-DLOG(DSQRT(PTK*ATP*ETA/K)/(KTP/K)**2) 
SIG=L53D-2*(KTP/K)**(3./2.)/COULOG 
U=(I*E*G/(2*PI*SIG*R*J))*a  +  Z)/a+(E*G/J)*(L+Z))- 
1      BETA*(KTP/E)/(1.+(E*G/J)*(1.  +  Z)) 
VT=VK+U 
WRITE(6,910)   TK,KTP/E,VK,a,VT 
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WRITE(6,920)   PTK/1.01D5, ATP, ETA, JI/J.R 

WRITE(6,930)   F.JE.JI.J 

WRITE(6',940)   COULOG,    SIG,    G 

WRITE(7,910)   TK,KTP/E,VK,U,VT 

WRITE(7,920)   PTK/1.01D5, ATP, ETA, JI/J.R 

WRITE(7,930)   F,JE,JI,J 

WRITE(7,940)   COULOG, SIG, G 

TK=TK+DTK 
20  CONTINUE 

STOP 
900  FORMAT(1X,I3,5X,F4.0) 

910  FORMAT(1X,F5.0,5X,F6.2,5X,F6.2,5X,F6.2,5X,F6.2) 
920  FORMAT(1X,F4.1,5X,F6.4,5X,F6.4,5X,F6.4,5X,1P,D10.4 
930  FORMAT(1X,1P,D10.4,5X,D10.4,5X,D10.4,5X,D10.4) 
940  FORMAT(1X,F4.1,5X,F7.1,5X1P,D10.4) 

END 
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